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FOREWORD 


SAAM  Is  a digital  computer  program  developed  for  the  analysis  of 
data  In  terms  of  models.  It  permits  simulation  euxd  data  fitting,  and 
contains  various  technlq]ue8  encountered  In  model  building. 

Although  developed  primarily  for  biological  systems  emd  more  spe- 
cifically for  kinetic  models,  the  program  Is  of  general  utility.  It 
differs  from  other  simulation  sad  analysis  systems  In  that  the  "language" 
Is  geared  towards  the  bio-medical  "system"  Investigator  and  Its  elements 
are  direct  counterparts  of  techniques  and  conceptualizations  xised  by  the 
experimenter. 

Model  building  Is  complicated  and  requires— In  addition  to  Intuition 
and  speculation — knowledge  of  mathematical  and  statistical  procedures  and 
their  limitations.  This  manual  Is  only  a brief  description  of  the  proce- 
I dures  used  In  SAAM  and  some  of  their  limitations.  For  additional  back- 
ground material  the  reader  Is  referred  to  the  reference  section. 

SAAM  Is  a large,  complex  program  and  Is  continuously  being  extended 
and  revised.  like  any  large  program,  SAAM  Is  difficult  to  coaqpletely 
debug,  and  It  probably  contains  some  undetected  errors,  even  through  It 
has  been  in  use  since  1939*  It  Is  recommended,  therefore,  that  the  user 
run  some  test  problems  of  his  own,  the  answers  to  which  he  knows.  We 
also  Invite  users  to  call  to  our  attention  any  questionable  results  which 
may  be  attributed  to  the  xn:ogram  and  not  to  errors  In  the  data. 

This  manual  Is  for  the  SAAM  23  version  of  the  program.  Revisions  and 
updates  will  appear  occasionally,  euid  will  be  sent  to  those  who  request 
that  their  names  be  placed  on  our  mailing  list.  A new  version  of  SAAM  Is 


novr  under  development  euid  vlll  be  known  as  SAAM  2k,  SAAM  2k  vlU  contain 
major  revisions  and  many  new  features,  and  a ccmpletely  nev  manual  vlll 
be  written  for  it. 
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INTRODUCTION 

SAAM  is  a general  purpose  computer  program  designed  to  fit  physical 
or  mathematical  models  to  data  hy  adjusting  the  parameter  values  of  the 
model  until  a '*best"  fit  is  obtained.  Any  set  of  mathematical  equations 
(differential,  integral  or  algebraic)  or  functions  may  serve  as  a model 
provided  an  ajialytical  or  numerical  procedure  exists  for  its  solution. 

An  open-ended  library  of  model  types  is  incorporated  within  the  program 
for  routine  use.  A partial  list  of  the  model  types  included  in  the  SAAM 
library  is  given  in  Section  II. 

The  program  uses  a ccanmon  data  input  format  for  all  types  of  models. 
This  is  made  possible  through  the  use  of  a single  set  of  operational  ele- 
ments in  the  program  and  a defined  equivalence  between  these  elements  and 
the  elements  of  each  model  type.  SAAM  permits  the  acceptance  of  "raw” 
experimental  data.  The  number  of  entries  required  for  the  specification 
of  models  and  constraints  has  been  minimized  to  simplify  use. 

As  will  be  discussed  more  fully  later,  SAAM  contains  a number  of 
features  designed  to  aid  an  investigator  in  his  model  building  efforts. 

As  a result,  the  program  has  become  quite  large. 

The  SAAM  22  version  of  the  program  is  in  FORTRAN  II  (with  a few  sub- 
routines in  FAP),  contains  over  200  subroutines  and  about  15,000  FORTRAN 
statements.  It  is  compiled  under  the  IBM  independent  FORTRAN  monitor 
for  routine  execution  in  an  IBM  T09*<-  using  the  CHAIN  facility.  A CHAIN 
tai)e  generated  from  the  relocatable  binary  links  is  saved  for  routine  use. 
At  execution  time  the  preaddressed  CHAIN  tape  is  mounted  and  a short  MAIN 
program  is  loaded  to  call  in  the  first  link.  Subsequently,  particular 
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links  are  brought  into  core  as  each  problem  requires. 

The  SAAM  22  version  has  also  been  converted  to  FORTRAN  IV  and  com- 
piled for  routine  use  on  the  CDC  56OO  and  CDC  66OO  systems. 

The  SAAM  23  version  is  a revised  version  of  SAAM  written  in 
FORTRAN  IV.  Most  of  the  modifications  are  minor  ajid  need  not  concern 
the  user.  Problem  decks  run  with  SAAM  22  will  also  run  with  SAAM  23. 

Problems  to  be  run  are  loaded  behind  the  MAIN  program.  Further 
details  on  the  orgemization  of  SAAM  into  blocks  are  presented  in 
Section  VIII  of  this  manual. 
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COMPUTATIONAL  STRUCTURE 

Formally  we  can  describe  a system  or  a model  by  its  responses,  a 


set  : 


t) 


where  F^  is  the  functional  relation  of  f^  to  a set  of  parameters  . 

X(Xi,  etc,),  some  set  of  condition  values  f^(fQ^^,  f^g^  •••) 

that  characterize  constraints  (such  as  initial  or  boundary  conditions), 

and  seme  independent  variables  t(ti,  tg,  .,,). 

Experimental  data  (q.°)  are  observables  of  the  system  and  as  such  are 

J 

estimates  of  some  theoretical  values  q.  that  are  functions  of  the  f. : 
1 

In  these  terms,  the  fitting  of  a model  to  data  implies  the  derivation  of 

a set  of  parameter  values  for  which  the  q will  ’^est"  fit  the  data,  q°. 

J J 

This  involves  three  stages  of  computations : 

1)  equations  solution  - to  calculate  f^  given  a set  of  values  for 

the  X and  f . 

o 

2)  matching  - to  convert  the  f^  to  the  theoretical  q^, 

5)  parametric  fitting  - to  adjust  the  X values  until  the  q^  '^est" 
fit  the  data  q?. 

When  the  functional  relation  F^  and/or  the  number  of  parameters,  X, 
are  \mknown,  the  fitting  of  a model  to  data  also  implies  derivation  of 
an  F^  and  a set  of  X*s,  This  is  referred  to  as  model  building,  and  will 
be  considered  only  briefly  in  this  manual, 

1,  Solution  of  equations 

This  stage  calculates  f^,  given  a set  of  X and  f^. 


The 
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computations  depend  on  the  type  of  model  and  on  the  availahle  ccmxnita- 
tional  methods.  For  example,  differential  equations  are  solved 

numerically  using  a 4 order  Runge-Kutta  method  (7)*  Linear 
differential  equations  vith  constant  coefficients  can  also  be  solved 
using  an  eigenvalue -eigenvector  method  or — in  special  cases--using  an 
analytic  method.  Each  method  of  solution  is  executed  by  a sep- 

arate subroutine.  Each  subroutine  calculates  values  for  f^  required  for 
the  calculations  of  any  q^  requested  in  a problem.  The  q^  requested  may 
correspond  to  observations  q?  or  to  any  simulated  quantity  of  interest. 

The  subroutines  for  the  solution  of  different  model  types  are  stored 
in  a library  of  ''models".  During  the  execution  of  a particular  problem 
the  appropriate  subroutines  are  brou^t  into  core  memory. 

2.  Matching 

An  observation,  q°,  is  a meas^ure  of  some  function,  q.,  of  the 
J J 

set  of  solution  values  f^  that  characterize  the  system  or  model.  Due  to 

experimental  constraints  the  q.  need  not  correspond  directly  to  the  f „ . 

J t 

Althou^  various  functional  relations  between  the  q^  and  f^  are  possible, 
at  present  only  linear  relations  are  provided  for: 


I 


"ji^i 


The  f^  required  for  the  calculation  of  the  q^  are  computed  and 
stored  in  the  equations  solution  stage.  If  the  are  known  from  the 
experimental  set  up,  the  q.  can  be  calculated  directly.  If  the  are 
not  known^  their  values  nay  be  estlnated  from  the  observed  <iuantltles 


i 
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q?  = K.\  CT..f. 
J J Z_.  1 


using  a least  squares  linear  regression  analysis  (8),  provided  the  num- 
ber of  independent  observations  is  equal  to  or  greater  than  the  nvimber  of 

unknovm  a...  With  the  a.,  thus  determined,  the  q.  can  be  calculated. 

Ji  Ji  J 

In  the  structure  of  SAAM  the  a.,  appear  as  secondary  parameters, 

as  distinguished  from  the  which  are  primary  parameters. 

The  calculation  of  a specified  set  q.,  given  a model  and  values 

J 

for  X,  f^,  amd  a^  is  referred  to  as  simulation  ajid  is  equivalent  to  the 

same  process  performed  on  an  analog  computer. 

5.  Rirametric  fitting  of  data 

This  stage  involves  the  adj\istment  of  the  model  parameters, 

until  the  calculated  q.  "best"  approximate  the  observations  q?.  In 

J J 

general,  the  q.  are  non-linear  functions  of  the  X,  and  a non-linear 
J 

least  squares  fitting  procedvire  is  employed.  Starting  with  initial 
estimates  for  the  parameter  values,  ope  calculates  a set  of  q^.  Adjust- 
ment of  the  initial  X values  is  then  mauie  using  a first  order  approxima- 
tion. 


"j 


Z-«i 


The  partial  derivative  of  each  q,  vith  respect  to  \ is  obtained  as  an 

J 

approximation  by  calculating  a Aq,  for  a small  change  AX.  and  setting 

J ^ 

8X, 


The  Sq^  are  approximated  by  the  difference  between  the  observed  and 
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calculated  q 


J 


= 


so  that 


I 


This  constitutes  a set  of  equations  from  which,  using  lineeur  regression 
analysis,  a set  of  6X^  can  be  calculated  to  minimize  (q°  - Since 

q^  are  non-linear ly  related  to  the  the  above  procedure  is  only  an 
approximaticxi  to  a least  squares  fit  and  is  iterated  until  convergence 
is  achieved. 

In  the  neighborhood  of  a least  squ£u:es  fit  this  procedure  yields  a 
variemce -covariance  matrix  for  the  X.  From  this  matrix  are  calculated 
estimates  for  the  standard  errors  of  the  X and  their  correlation  coef- 
ficients. These  estimates,  however,  do  not  necessarily  reflect  the  true 

xuacertainties  in  the  X*s  due  to  non-linearities  in  the  q.  with  respect 

J 

to  the  X^. 
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MODEL  TYPES  IN  PROGRAM  LIBRARY 

Different  types  of  models  cam  be  processed  by  the  program.  Each 
type  of  model  is  coded  so  that  the  data  deck  will  be  correctly  interpreted. 
To  simplify  the  preparation  of  data  a common  nomenclature  and  data  form 
have  been  adopted  for  all  types  of  models,  amd  equivalences  are  defined 
between  the  parameters  of  the  model  type  and  the  common  ncmienclature . 
Operational  units  and  graphical  symbols  are  Introduced  as  aids  in  describ- 
ing the  equivalences.  A detailed  discussion  of  each  model  type  is  given 
in  Section  VI  on  Model  Types. 


II-6 


May  1966 


Program 

Element 

Description 

Symbol 

f^Ct) 

function  ^ 

© 

f^io) 

initial  value  of  f^(t)  at  t = 0 

canposite  function  k,  a function  of  one 

or  more  f . *s 
1 

A 

primary  parameter 

secondary  parameter:  summing  coefficient 

0^^ 


K. 

1 


secondary  parameter:  proportionality  coefficient 

% = 


K.f. 

1 1 


or 

= ''i  Z 

F.  fuD,ction  generator  j 

J 

t independent  variable 

9 second  independent  variable 

matrix  i 
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j INPUT  FORMAT 

! A number  of  headings  are  used  in  setting  up  a problem.  These 

j are  described  in  this  section  in  the  order  required  by  the  program, 
; The  formatting  of  information  under  each  heeuiing  makes  use  of  two 
I types  of  fields:  integer,  and  decimal. 

I Integer  fields  require  that  no  decimal  point  be  entered,  and 

I that  the  number  be  at  the  right  end  of  the  field. 

I 

Decimal  fields  permit  the  number  to  be  entered  anywhere  in  the 
indicated  field  but  the  decimal  point  must  be  entered,  even  for  a 
whole  number. 
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The  problem  deck  is  organized  in  sections  which  must  be  in  the 
following  order: 

A.  Required  sections: 

1.  card  1 - clock  option  and  problem  I.D. 

card  2 - nximber  of  components  and  number  of  iterations 
card  3 ~ timing  and  convergence  criteria 
card  k - options  and  model  type 

2.  data 

*26' (termination  card*) 

3.  initial  conditions 

*26* 


4.  Kappas 

*26* 


3 . lambdas 

*26* 


6.  Sigmas 

*26* 


7.  Dependence  Relations 

•26 » 


8,  Statistical  constraints 

*26* 


B.  Additional  sections  as  required  by  the  problem: 


1.  T-interrupt  changes  in  f^ 

•26  • 

•26* 

2,  T-interrupt  changes  in  parameters. 

•26 » 


3.  Normal  equations  entries 


C.  Special  entries 

These  are  required  only  for  particular  model  types.  How 
and  where  they  are  to  be  entered  are  described  under  model 
type. 

♦Footnote  - for  required  sections  (2  throu^  8)  the  termination  cards 

*26*  must  be  present  whether  or  not  other  cards  are  entered. 
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LIMITS  IMPOSED  BY  PROGRAM 

Because  of  physical  limitations  in  the  computer  the  following 


at 

present,  imposed  by  the  SAAM  program: 

of 

components 

< 

25 

of 

primary  and  secondaiy  parameters 

66 

of 

dependence  relations 

< 

of 

variable  paarameters 

25 

of 

data  and  "statistical  constraints" 

250 

of 

T-interrupts 

4 

CARD  1 


o 

CO 


tr\ 

C— 


(M 


OJ 

CU 


io» 

rH 


OJ 

r— 1 
+ 
Q 


00 
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CARD  1 


Coliomn  1.  Enter  a "l"  or  "2"  to  clue  the  prograim  that  this  is  the 
beginning  of  a problem  deck. 

The  "2"  is  used  when  an  ON-LINE  clock  is  available  on  the 
ccmiputer.  This  will  result  in  a printout  of  the  starting  time 
for  a problem  and  its  elapsed  execution  time.  A "l"  in  this 
field  will  bypass  the  clock. 

Colxomn  5 ->  8.  Enter  '^AAM  25”  to  identify  the  version  of  the  program 
used. 

Column  10  •>  21.  Problem  identification  number:  three  initials 

followed  by  up  to  nine  characters  (numbers,  letters,  decimals). 

Column  25  72.  Users  name  and  description  of  problem. 

User  is  free  to  enter  whatever  he  wishes  in  this  field. 

It  will  be  printed  in  the  output  as  is. 


Card  2 


o 

CO 

f 

tOy 

O 

+ 

VO 

VO 

trv 

vO 

cH 

vO 

O 

vO 

+ 

vO 

LTV 

LT\ 

iT\ 

f 

rH 

LTN 

o 

LTN 

ON 

-=i- 

to 

«t-t  O 
O >H  <U 
-P  tao 

o 

to 

• cj  0) 
O U ^ 

a (0  0 

■P  *H 
H 

ov 

cu 

No.  of 
Components 
(integer) 

o 

<M 

ON 
1 — 1 

o 

fH 

VO 

f 

CVJ 

iH 
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CARD  2 


Column  19  ■*  20.  Number  of  components.  This  corresponds  to  the 

largest  index  associated  with  any  component  of  the  model. 

Column  29  -*■  50*  Maximum  number  of  iterations  for  least  squares 

convergence.  Computation  may  be  terminated  internally  before  this 

maximum  is  reached  if  convergence  criteria  are  satisfied. 

"Blank"  or  "zero"  entry  in  this  field  is  interpreted  as 

ZERO  ITERATIONS  and  results  in  just  one  solution  based  on  initial 

values  of  primary  parameters.  With  secondary  parameters  fixed 
is 

this  equivalent  to  a SIMULATION. 


Card  3 


s 

t 

1 1 

CONMIN 

(decimal) 

o 

t 

iH 

VO 

o 

t 

rH 

Lf\ 

E 

(decimal) 

0 

ITS 

1 

P 

1 

(decimal) 

s 

t 

•H 

rOy 

o 

% 

P 

(decimal) 

f 

fH 

CVJ 

11— > 20 

f 

TIME 

FACTOR 

(decimal) 

2 — »10 

rH 
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CARD  3 


Column  2 -►  10.  "Blank”  unless  the  internally  allotted  time  for  a 

single  solution  of  the  differential  equations  is  to  be  modified. 

Internally  allotted  time  is  O.5  units.  (Computational  units) 
Enter  smaller  or  larger  value  as  desired. 

Column  21  -»■  30*  "Blank".  See  Section  VII-5  for  eidditional  information. 
Column  31  Enter  ”.01".  This  entry  indicates  fractional  change 

in  parameter  X used  to  calculage  dq/^X.  Usual  entry  .01. 

Column  kl  -*•  51»  This  is  convergence  test  value.  Usual  entry  .98. 

Column  61  -►  70.  This  is  also  convergence  test  value.  Usual  entry  .98» 


Card  U 


o 

00 

+ 

KN 

O 

C^ 

+ 

MD 

MD 

LPl 

+ 

iH 

MD 

O 

VD 

+ 

VO 

Lf\ 

LT\ 

LTN 

i 

r— I 

ir\ 

J 1 

w w (U^ 
n -p  ^ 

O 

Ll^ 

Q o C 0) 
S O M tjO 

ON 

INTER- 

MEDIATE 

RESULTS 

VO 

A 

MATRIX 

LTN 

PLOT 

-d- 

PARTIALS 

MATRIX 

KV 

COVARIANCE 

MATRIX 

C\J 

rH 

Column  2. 

"Blank*'  - no  print  or  punch  of  covariance  matrix, 
"1"  - print  and  punch  covariance  matrix. 

Column  3. 

"Blank"  - ^ printout  of  partials  matrix, 
"l"  “ printout  of  partials  matrix. 

Column  1+, 

Plotting  options : 

"Blank"  - no  plots 
"1"  - semi “10^5  2 page  plot, 

"2"  - semi -log  1 page  plot, 

"5"  - linear  2 page  plot, 

"i|."  - linear  1 page  plot. 

Colxjmn  5« 

"Blank"  - ^ print  or  punch  of  A matrix, 
"l"  - print  and  punch  A matrix. 

Column  6. 

"Blank"  - no  print  of  intermediate  results, 
"l"  - print  intermediate  results. 

Column  k9 

► 50*  Type  of  model  solution, 

"Blank"  - program  will  choose  automatically 
among  solutions  1,  2 and  4. 
integer  entry  - nxomher  corresponding  to  desired 
model  type.  (See  code  in  section  on  model  types  #) 

DATA 
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DATA 


This  heading  permits  the  specification  of  functions  to  be  calcu- 
lated for  simulation  or  for  the  fitting  of  data.  The  function  numbers 
(i),  the  independent  variables  (t,9),  observed  values  (q^)  and  the 
statistical  wei^ts  are  entered  as  follows ; 

colvunn  2 -*■  Conponent  number  - specifies  subscript  i of 

(function  or  summer)  as  defined  by  the  model  code, 
coltunn  15  25.  Value  of  independent  variable  t, 

column  27  **-  40,  Observed  value  q^(t,9),  if  available, 
column  k2  55,  Relative  statistical  weight  of  observed  value, 
column  57  Blank 

column  59  ♦ 72.  Independent  variable  9. 
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Each  line  lander  this  format  constitutes  a request  for  a single 
computation  q^(t,9).  The  execution  sequence  of  computational  requests 
may  he  changed  internally  hy  the  program. 


Control  Cards 


1)  Entry  control  card  (EC) : modifies,  defines  or  interprets  the 

entiles  on  the  data  cards  that  follow  it.  Control  is  terminated 


by  a new  entry  control  card  or  by  the  termination  of  entries  ~\ 
xmder  the  DATA  heading.  A control  card  is  not  considered  a datum,' 

I 

column  2 -*■  5*  Entry  of  100.  + X.  Designates  an  EC  card  | 

I 

and  assigns  X to  all  succeeding  entries  under  COMP,  \ 

i 

except  when  X = 0.  When  X = 0.  COMP  entries  remain  ! 

xinchanged.  For  example,  an  entry  IO5.  will  assign  a j 

CCMP  No.  5*  to  all  succeeding  entries;  a 100.  will  leave 
them  unchanged. 

column  13  -*■  25.  An  entry  X in  this  field  will  add  X to  all 
entries  that  follow  in  this  field.  This  provides  for  a 
shift  in  the  "t"  coordinate.  X may  be  positive  or 
negative . 

column  -*•  kO,  An  entry  X in  this  field  multiplies  by  X 
all  entries  that  follow  in  this  field.  When  X = 0,  or 
blank,  the  entries  that  follow  remain  unchanged, 
column  14-2  -►  55  and  column  57  are  used  jointly  as  a code 
for  statistical  weight  assignment  as  follows : 


Pe 


1 1, 
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Entry  Control  Card 

Interpretation  of  weight  field  (columns  42  -►  55) 
of  succeeding  cards 

Columns 
42  --  55 

column 

. 51 

blank 

blank 

relative  wei^t 

blank 

1 

standard  deviation  of  observations 

X 

1 

standard  deviation  of  observed  value  equals  X. 
data  card  wei^t  field  Ignored 

X 

2 

coefficient  of  variance  of  observed  value  equals  X. 
data  card  weight  field  ignored 
X must  be  non  zero 

Note : 1)  Weights  are  calculated  and  normalized  by  the  program  vhen 

I 

the  standBurd  deviations  or  coefficient  of  variance  are  given. 


2)  SAAM  assigns  zero  wei^t  whenever  the  given  or  calculated 
standard  deviation  is  zero. 

3)  The  weight  assignments  for  a given  problem  must  be  compatible. 
Either  xelative  weights  or  wei^ts  betsed  on  standard  deviations  must  be 
used  exclusively.  Both  types  of  wei^tlng  will  not  be  accepted  in  a 
single  problem. 

j 4)  See  section  on  Methods  for  discussion  of  statistical 

! wei^ts  and  the  calculation  of  SIG. 

1 

2)  t -interrupt  control  card  (TC) ; interrupts  the  computations  euad 
permits  changes  to  be  introduced  during  the  solution  of  a problou. 
Specifically,  at  any  specified  "t",  existing  parameter  values  may  be 

! switched  to  new  values,  and  Af^  may  replace  or  be  added  to  the  values 

1 

j f^  existing  at  '*t”. 


I 


I 
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colunm  2 -*■  5»  'Hie  entry  ”126,”  signifies  a TC  card 

effective  at  the  t-value  assigned  to  the  datum 
immediately  preceding  it.  This  remain  in  effect  until 
another  TC  card  or  termination  card  is  encountered, 

(EC  cards  work  independently  of  the  TC  cards,) 
column  12  25. 

a)  A blank  in  this  field  Indicates  that  the  t-values 
of  the  data  that  follow  are  to  be  interpreted  with 
respect  to  a new  t scale  starting  with  t = 0 at  the 
TC  card. 

b)  A ”1,”  in  this  field  indicates  that  the  t-values  of 

the  data  that  follow  refer  to  the  same  t scale  as 
those  preceding  the  TC  card.  ^ 

column  27  40, 

a)  A blank  in  this  field  indicates  that  a new  set  of  f, 
values  replace  the  current  fj^  values  at  the  t-interrupt. 
The  new  values  are  specified  under  ”t -interrupt  changes 
in  f^.” 

b)  A ”1,”  in  this  field  Indicates  that  a set  of  values 
Afj^  (entered  tuiAer  ”t -interrupt  changes  in  fj^”)  are 
to  be  added  to  the  existing  f^^  at  the  t-interrupt. 

colvuan  42  -*■  55*  Blank, 
column  ^7*  Blank 
column  59  72.  Blank 

5)  Data  Generation  Control  Card,  (GC) : is  used  as  an  expedient 

to  generate  ”axtificial"  data.  This  control  card  generates  X 


< 
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data  entries  at  intervals  At  staurting  vith  the  t-value  of  ^he 
last  datum  (entered,  or  generated).  The  component  number  and 
velght  of  the  last  datum  (entered  or  generated)  is  carried 
over  to  the  generated  entries.  A ntjmher  of  GC  cards  can 
follovr  each  other. 

column  2 5«  Tlie  entry  ”200.'*  signifies  a GC  ceurd 

column  12  -*■  25.  Enter  At,  the  t-intesrval  between  generated 
data  points, 
column  27  ■*  Blank 

col\imn  42  55*  Enter  X,  the  number  of  data  points  to 

be  generated, 
column  37*  Blank 


column  59  •*"  72.  Blank 


INITIAL  CONDITIONS 
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This  heading  permits  the  entry  of  f^(0)  values  for  each  component  i 
Only  non-zero  entries  need  he  made. 


column  4 -»■  5 

component  number,  i 

column  12  25 

initial  conditions  f^(0)  ' ‘ 

colixmn  42  ->  55 

V 

column  56  ->  TO 

U 
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PARAMETERS 

SAAM  recognizes  three  types  of  parameters:  KAPPAS,  LAMBDAS  and 

SIGMAS.  The  lAMBDAS  are  usually  non-linear ly  related  to  the  functions 
f^(t)  and  req.uire  iterative  adjustment  for  a least  squares  fit.  The 
SIGMAS  and  KAPPAS  are  linear  functions  of  the  observations.  The  hierarchy 

- i 

of  canpu tat ions  in  the  program  is  as  follows ; 


All  parameters  are  doubly  subscripted,  although  the  second  subscript 
of  Kappas  (always  zero)  is  frequently  omitted. 

From  the  point  of  view  of  adjustment  of  parameters  to  fit  data, 
the  parameters  are  classified  as: 


assigned  in  advance. 

dependent:  parameter  takes  on  a value  as  prescribed  by  a 

dependence  relation. 

adjustable : parameter  is  free  to  adjust  as  required  by  the 

program  in  data  fitting. 

Fixed  parameters  may  change  during  a solution  but  always  in  a pre- 
determined manner,  independent  of  the  fit.  Such  changes  may  occur,  for 
example,  in  connection  with  a t-interrupt. 

Dependent  parameters  are  usually  dependent  on  other  parameters  and/or 
constamts  and 

are  indirectly  adjusted  in  connection  with  the  fitting  of  data, 

A parameter,  once  classified  as  dependent,  remains  so  throu^out  the 


fixed : 


the  parameter  takes  on  preset  values 
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solution  emd  always  obeys  the  same  functional  dependence  relation 


A parameter  is  designated  as  dependent  by  a code  on  its  entry 
card  (“l"  in  column  60).  The  dependence  relation  is  specified  under 
a separate  entry  (Dependence  relations). 

At  present  dependence  relations  are  restricted  to  linear  relations 
only,  such  as 


the  right  side  of  the  equation  can  be  fixed,  dependent  or  adj\istable. 

The  program  resolves  these  dependence  relations  internally,  so  that 
dependent  parameters  are  finally  expressed  in  terms  of  adjustable  and 
fixed  parameters  only. 

Whereas  sigmas  and  kappas  can  be  dependent  on  other  sigmas,  kappas, 
lambdas,  or  constants,  lambdas  can  be  dependent  only  on  lambdas  and 
constants . 

Adjustable  parameters  have  upper  and  lower  limits  assigned  to 
them.  These  limits  are  never  violated  in  adjusting  the  parameter  during 
data  fitting.  The  program  identifies  a parameter  as  adjustable  because 
its  upper  limit  is  greater  than  its  lower  one. 

An  adjustable  parameter  is  entered  with  an  initial  estimate*  emd 
a maximum  and  minimum  limit.  In  addition,  an  estimate  of  a standard 
deviation  can  also  be  entered,  when  known  a priori.  Such  information 
may  be  known  frcxn  sources  other  than  the  data  to  be  fitted.  The  standard 

♦Initial  estimates  are  not  required  for  sigmas  and  kappas  unless  a 


k,h 


The  parameter  on  the  left  side  is  the  dependent  one  and  the  on 


standard  deviation  is  given  for  them 
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deviation  is  associated  vith  the  initial  estimate  value  and  is  combined 
statistically  with  the  data  by  the  program  to  derive  best  estimates  for 
the  adjustable  parameters.  More  extensive  statistical  constraints  can 
also  be  entered  vmder  the  separate  heading  ’’STATISTICAL  CONSTRAINTS”. 

It  may  be  noted  that  the  statistical  constraints  on  a parameter 
are  independent  of  the  upper  and  lower  limits.  The  latter  are  not 
involved  in  calculating  corrections  for  the  adjustable  parameters  but 
are  used  to  limit  the  corrections  to  within  the  specified  region. 

Parameters  can  also  be  f -dependent . 'IMs  means  that  during  compu- 
tations their  values  are  modified  by  an  f^(t)  for  every  t.  f -dependence 
is  restricted  to  two  forms : the  parameter  value  is  multiplied  or 

divided  by  some  f^(t).  At  present  only  model  type  4 permits  f -dependence . 

The  kappas  are  read  first  by  the  program.  Before  reading  this 
category  the  program  automatically  assigns  fixed  values  = 1.  for  all 
components  of  the  model.  These  values  are  overwritten  by  the  entries. 
Therefore  only  / 1 need  be  entered. 

The  lambdas  are  entered  next  followed  by  the  sigmas.  Only  lambdas 
and  sigmas  required  in  the  model  need  be  entered. 

The  order  of  the  entries  tmder  each  heading  is  not  critical. 
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column  4 -►  5»  1 subscript  of  parameter. 

colmnn  9 10,  2™"  subscript  of  parameter.  For  KAPPAS  second 

subscript  is  always  zero. 

column  12  -»•  25.  Value  of  parameter.  This  field  is  ignored  by  the 
program  for  variable  SIGMAS  and  KAPPAS  unless  STANDARD  DEVI- 
ATION field  (62  -*-72)  has  a,  non-zero  entry, 
column  27  •>  40,  Lower  limit  for  euijustable  parameter, 

column  42  - 55*  Upper  limit  for  adjustable  parameter.  Upper  limit 

must  always  be  greater  than  lower  limit  for  all  adjustable 
parameters . 

column  57  59 • Blank;  if  there  is  no  f -dependence . 

+ i:  if  parameter  is  to  be  multiplied  by  f^(t) 

throu^out  the  computation, 

- i:  if  parameter  is  to  be  divided  by  f^(t) 

throu^out  the  computation. 
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coliunn  60  : Blank:  if  parameter  is  fixed  or  adjustable. 

”1”:  if  parameter  is  dependent.  Dependence 
relations  are  entered  separately  under  '*DEEENDENCE  RELATIONS '* 
column  62  -►72:  enter  X for  adjustable  parameter  only  when  the 

value,  P,  of  the  parameter  is  known  with  a standard  deviation 
+ X : value  = P + X. 


Y OF  INFORMTION  TO  BE  ENTERED  ON  A CARD  DEFINING  A PARAMETER 


DEPENDENCE  RBIATIONS 
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DEIENDENCE  RELATIONS 


Dependence  relations  are  limited  to  linear  dependences  of  the  type 

+ Aqq  for  specified  k and  t 

where  the  Ej^  are  fixed,  dependent  or  adjustable  and  with  the  further 

constraint  that  \*s  can  be  dependent  only  on  X*s. 

column  k -►  5»  Enter  ^ for  dependent  ^ or  k.  . 

1 J 1 J 3- 

coliunn  9 10*  a)  Enter  j for  dependent  cr^^. 

b)  Leave  this  field  blank  for  K. . 

coliamn  19  20, 

a)  Blank  for  dependence  on  a constant, 

b)  Enter  i for  dependence  on  X^^  or  K^, 
column  2k  -»•  25 

a)  Blank  for  dependence  on  a constant  or  K^, 

b)  Enter  ^ for  dependence  on  Xj|^j  or 
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column  27  ■*  ^0 

a)  Enter  the  coefficient  of  the  depended-on  i>araineter 

b)  Enter  a depended-on  constant 
For  example: 

= .5  *"  *^15  ^1^0 

Each  term  of  the  dependence  relation  is  entered  on  a separate 
card  as  follows : 


k 

? 

9 

10 

19 

20 

2k 

25 

27—*  40 

2 

5 

6 

T 

.5 

2 

3 

1 

-.7 

2 

5 

1 

0 

16. 

2 

5 

0 

0 

12. 

The  subscript  (0,0)  is  reserved  for  a constant. 
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STATISTICAL  CONSTRAINTS 


STATISTICAL  CONSTRAINTS 


Information  is  entered  under  this  heading  when  statistical  uncer- 
tainty associated  with  a parameter  (lambda,  sigma  or  kappa)  or  linear 
COTbination  of  parameters  is  known  independently  of  the  data.  This  is 
a more  general  procedure  than  entering  standard  deviations  directly  in 
columns  62  -►  72  of  the  parameter  cards.  The  latter  procedure  is  used 
merely  for  convenience  when  the  statistical  constraints  apply  to  a 


single  parameter  only  and  when  its  estimated  value  agrees  with  the  initial 


I estimated  value  and  X a standard  error  of  the  estimate,  the  general 

I 

form  for  statistical  constraint  may  be  written  as 


This  is  entered  on  a number  of  cards,  one  for  each  term.  Y + X is 


value  employed  in  the  fitting. 


) 


Y + X 


(for  i,j  involved) 


entered  with  the  last  term, 
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colvmin  U -»•  5.  Enter  "i"  for  or 

column  9 -*■  10, 

a)  Enter  "j"  for  X.  . or  a.  .. 

XJ 

b)  Blank  for 

column  12  -»■  25.  The  coefficient  (A^j ) of  the  parameter  entered  in 
colxmm  4 -*■  10, 

column  26  -»•  40,  The  constant  (Y,)  on  card  of  last  term  only, 
column  4l  55*  The  error  (X,)  on  card  of  last  term  only. 


For  example: 

•T  ^13,4  ^1  ” 5 i 


sa 

4 

5 

9 

10 

12  25 

27  -►  40 

42  55 

1 

2 

.7 

1 

4 

-.5 

1 

0 

1, 

5. 

1, 

4 

T” INTERRUPT  CHANGES  IN  f. 


f 


f 

f 

I 

S 
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T-IKTERRUPT  CHANGES  IN  f. 

1 


Changes  in  associated  with  t-interrupt  control  cards  are  entered 
under  this  heading.  The  amounts  entered  will  he  added  to,  or  he  used  in 
place  of,  the  existing  value  at  the  interrupt,  as  instructed  hy  the 
TC  caxd  (see  instructions  on  TC  cards).  Two  "26”  termination  for  this 
category  are  required  whenever  any  t-interrupt  (t)  cards  are  involved, 
regardless  of  whether  or  not  changes  in  f^  are  to  he  made, 

column  4 -*•  5»  The  number,  i,  of  the  component  whose  f^  is  to 

he  changed, 

column  7 -►20,  Amount  of  change  at  first  t-interrupt  (TCI) 

column  21  -►  35 • Amount  of  change  for  second  t-interrupt  (TC2) 

column  56  -►  50 • Amount  of  change  for  third  t-interrupt  (TC3) 

column  51  -►  65,  Amount  of  change  for  fourth  t-interrupt  (TC4) 

If  no  entries  are  made  under  this  heading  for  a particular  component 
^ zero  entries  are  automatically  assumed  internally.  This  means  that 
"zero"  will  he  either  added  to  or  replace  f^,  as  instructed  hy  the  entry 
on  the  TC  card, 

TWO  (2)  "26"  termination  cards  are  required  after  the  last  entry 
under  this  heading. 


rPT  CHANGES  IN  P; 


PP 

+ 


ENTER 

10  for  TC  0 

11  for  TC  1 

12  for  TC  2 
15  for  TC  3 
Ik  for  TC  k 

ON 

P_— — _ 

PROBLEM 

NUMBER 

VO 

tv 

+ 

KN 

(integer) 

§ 

u 

o 

ITN 

ON 

-d- 

o 

p 

IT\ 

c 

H 

-d- 

B 

(decimal) 

o 

-4* 

T 

t- 

<M 

A 

(decimal) 

ITN 

CVl 

f 

a 

V* 

(integer) 

^ . 

10 

ON 

o 

l/N 

a 

H 

-4* 
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T-INTERRUPT  CHANGES  IN  PARAMETERS 


The  values  of  fixed  and  adjustable  lambdas,  and  fixed  sigmas  sind 
kappas  can  be  changed  at  each  t -interrupt. 

To  facilitate  the  description  and  limitations  of  this  feature,  we 
introduce  the  following  notation: 

X . . = value  of  parameter  i,J  prior  to  the  start  of  a solution 


for  f^(t). 


,th 


X.  .[TC(K)]  = value  of  parameter  i,J  as  set  at  the  K t-interrupt. 

Changes  can  be  made  at  each  t-interrupt  in  accordance  with  the 
following  relation: 


X^^rTC(K)]  = A Xj^^  + B 

X.  . and  X.  . can  represent  any  fixed  or  eidjustable  lambdas,  or  fixed 

i#J  n,* 

sigmas  or  kappas  and  A and  B are  arbitrary  constemts. 

If  any  parameter  is  changed  at  a t-interrupt,  TC(K),  all  other 
parameters  automatically  revert  to  their  values  preceding  the  very 
first  change;  Xj^^: 


X^/TC(K)l  * 
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If  no  parameter  is  changed  at  a TC(k),  then  all  parameters  retain 
their  values  prior  to  the  TC(k) 

X^^[TC(k)1  h Xj^^[TC(K-1)1  f or  k > 1 

^.?[TC(1)]  = 

A "26"  card  is  required  as  teimination  of  parameter  changes  for 
each  t -interrupt,  regeurdless  of  whether  or  not  actual  parameter  changes 
are  to  he  made. 

colrann  14-  -»•  5»  Enter  ^ for  X^^[TC(k)1 

column  9 -*■  10.  Enter  for  X.  .[TC(k)1 

^ J 

xolumn  12  -»•  25.  Enter  constant  coefficient  A 

column  27  -►  40,  Enter  constant  B 

column  Ml-  -►  i4-5 . Enter  h for  • 

column  k-9  -*•  50.  Enter  ^ for 


Note : It  should  he  noted  that  dependent  parameters  have  their  values 

recalculated  after  each  change  of  parameter  values  at  a t-interrupt. 
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NORMAL  EQJATIONS  ENTRY 

Frequently  it  is  desired  to  add  to  the  normal  equations  generated 
in  the  solution  of  a problem  einother  set  of  normal  equations  that  repre- 
sents the  aggregate  of  some  previous  information,  (For  example,  know- 
ledge of  a population  of  which  the  single  study  is  a member).  Two  options 
are  provided.  The  first  option  permits  the  entry  of  the  matrix  of  normal 

! equations  A and  a vector  C such  that  the  combined  solution  will  yield 
P P 

(A..  |-A)x  = C + l-O. 

P P ^ 

A and  C are  the  nomal  equations  components  normally  generated  by  the 

model  solution,  Z and  Z are  variances  associated  with  the  data  and  the 

P 

added  normal  equations,  respectively. 

Under  option  2,  the  added  normal  equations  satisfy  the  combined 
solution 

(A  + |-  A )x  = 0 + |-  A (X  - X) 

P P 

where  is  a vector  of  reference  values  belonging  to  A^  and  C^,  and  X 
are  the  values  of  the  adjustable  parameters  in  the  problem  at  the  time 
normal  equations  axe  calculated. 

Option  1; 

1)  Enter  "1"  in  column  65  of  card  4 of  the  data  deck, 

2)  Behind  the  last  card  of  the  deck  add  the  following  in  the 
order  indicated 

a)  a caxd  with  the  value  of  Z^  in  columns  21  -»•  55, 


E format 
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b)  The  matrix  is  entered  in  accordance  with  matrix 
entry  format  described  later.  Only  one  half  of  the 
symmetric  matrix  need  be  entered. 

c)  The  vector  C is  entered  as  a one  row  matrix,  in 

P 

accordance  with  matrix  entry  format. 

Option  2: 

1)  Enter  "2"  in  column  65  of  card  4 of  the  data  deck. 

2)  Behind  the  last  card  of  the  deck  add  the  following  in  the 
order  indicated 

a)  a card  with  2^  in  columns  21  -►  35  (E  format). 

b)  The  matrix  A - in  accordance  with  the  matrix 

P 

entry  format  described  later. 

c)  Xp  is  entered  as  a single  row  matrix. 

Note:  It  is  essential,  of  course,  that  the  order  of  the  adjustable  X*s 

be  the  same  in  the  problem  and  in  the  added  matrix  of  normal 


equations 
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MATRIX  ENTRY  FORMAT 


The  following  format  is  used  in  the  program  whenever  a matrix 
is  read  in  for  special  purposes: 

Row  dimension  of  matrix  (integer) 

Column  dimension  of  matrix  (integer) 

Row  of  first  matrix  element  on  card  (integer) 
Coltunn  of  first  matrix  element  on  card  (integer) 

Values  of  consecutive  matrix  elements 
(row -wise)  in  E Format,  (Blank  is 
interpreted  as  zero). 


^ 3 

It. 

8 

9 

15 

-»■  14 

18 

- 19 

21 

ro 

45 

^7 

- 58 

60 

^ 71 

Termination  of  matrix  input  is  indicated  hy  a card  with 
"26"  in  column  k. 
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SECTION  IV 


OUTRJT 


IV  - 1 


August  1966 


OUTPUT 

Routine  printout  is  produced  for  every  problem  deck 
submitted  to  SAAM.  Special  printout  is  added,  where  needed, 
for  certain  model  types.  Optional  output  is  also  available 
at  the  users  request.  This  includes  punched  as  well  as 
printed  output. 

Routine  printout  includes : 

1.  A listing  of  the  problem  deck  with  each  card 
printed  as  reformatted  by  SAAM. 

2.  SAAMs  reorganization  of  the  information  in  the 
problem  deck. 

5.  Parameter  values  and  a table  of  the  initial 
solution  values.  (Zero  iteration) 

k.  Results  for  each  iteration  associated  with  the 
fitting  of  the  data. 

5.  Final  results  which  include  parameter  values  euid 
the  corresponding  solution  table  for  the  '*best  fit”  with 
estimated  standard  errors,  and  correlation  coefficients. 

Optional  printout  includes  items  requested  under  OPTIONS 
entries  (card  4). 

Routine  emd  optional  outputs  are  described  in  greater 
detail  under  each  of  the  printout  headings.  These  are  given 
belovf  in  alphabetic  order. 

Special  outputs  in  connection  with  particular  model  types 
8ire  described  separately  at  the  end  under  the  appropriate  model 
type  or  subroutine. 
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A MATRIX  BEFORE  INVERSION  - contains  four  partitions  as  shown  below. 
Partition  I contains  coefficients  A(l,j)  of  normal  equations  for  lambdas 
only,  A * RES  = CR,  Kirtition  IV  contains  the  matrix  for  normal  equations  I 
to  calculate  sigmas  and  kappas 


I 

II 

III 

IV 

Partition  II  contains  matrix  of  coefficients  and  bK/b\,  Partition 

III  contains  zeroes. 

Option:  Column  5^  card  4 


C(l)  - denotes  i 


th 


component . 


CONAB  - is  a scaling  factor  for  the  RES  vector  to  achieve  a 
lower  sum  of  squares  The  value  of  CONAB  indicates  the  factor  by 

which  the  calculated  correction  vector  for  the  adjustable  lambdas  is 
multiplied. 


CORRECTED  SUM  OF  SQUARES  W PREVIOUS  ITERATION.  - Whenever  statis- 
tical constraints  are  given,  their  statistical  weights  are  renormalized 
in  each  iteration.  To  compare  the  sums  of  sqviares  of  two  iterations — 
as  a test  of  convergence— the  sum  of  sqviares  of  a previous  iteration  is 
recomputed  with  weights  assigned  in  the  current  iteration. 


CORRECTIONS  FOR  ADJUSTABLE  LAMBDAS  - lists  the  RES  vector  as  it 
will  be  used  in  the  current  iteration.  The  order  is  that  of  the  adjust- 
able lambdas. 
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j 

’ CORRELATION  COEFFICIENTS  - the  matrix  of  correlation  coefficients 

for  lambdas,  sigmas  and  kappas  is  partitioned  as  follows : 


X*s 

0 

0 

CT*s  and  K*s 

The  adjustable  parameters  appear  within  each  partition  as  ordered  lander 
PARAMETERS  but  with  the  lambdas  first,  then  Sigmas,  then  Kappas, 


COVARIANCE  MATRIX  - covariances  for  lambdas,  sigmas  and  kappas. 
Printout  is  partitioned  as  follows: 


covariance 
matrix  of 
lambdas 

0 

0 

covariance 
matrix  of 
sigmas  and 
kappas 

Option:  column  2,  card  k 


CR  VECTOR  - vector  elements  CR(I)  of  normal  equations  A * RES  = CR 
Option:  ccdtunn  card  k 


D - the  internal  order  of  data,  is  given  as  the  first  column  of 
the  solution  table.  Data  are  ordered  by  t values  within  each  block 
delineated  by  t-interrupt  control  cards. 
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FINAL  VALUES  - full  printout  of  values  for  the  '*best  fit"  arrived 
at  in  the  run.  Included  is  a table  of  dependent  and  adjustable  peuram- 
eters  with  their  estimated  standard  and  fractional  deviations,  and  the 
matrix  of  correlation  coefficients. 


INFORMATION  CONNECTED  WITH  CALCUIATION  AND  MODIFICATION  OF  RES  - 
the  RES  vector  as  solved  for  from  the  normal  equations  is  modified  and 
values  connected  with  this  modification  are  printed  in  table  form, 
(intended  as  an  aid  in  the  program  development). 


INITIAL  CONDITIONS  - is  a table  that  gives,  for  each  compcment 
(C(l)),  the  F(I,0),  the  type  of  component  (jUDY(l)),  the  changes  in 
F(l,t)  at  the  t-interrupt  (QP(i)),  smd  the  vectors  V(l)  and  U(l)  as 
specified  by  the  problem  deck. 


JUDY(i)  - is  a code  that  defines  the  type  of  component. 

JUDY(I)  = 1 means  that  component  I is  a function  f ^ (t ) . ^e  code  "2" 
signifies  a sxjjmning  component. 


K - represents  Kappa 


LITTLE  A BEFORE  MCOIFICATIC®  - is  the  triangular  matrix  a derived 
from  the  normal  equations  matrix  A: 

T 

a a = A 


Option:  column  5>  card  k 
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I 

j LITTLE  A- INVERSE  BEFORE  MODIFICATION  - is  the  inverse  of 

I "little  a",  derived  from  the  normal  equations  matrix  A 


Option:  colvimn  2,  card  4. 


PARTIAIS  OF  DATA  POINTS  WITH  RESPECT  TO  ADJUSTABLE  LAMBDAS  - lists, 
in  the  first  two  colvunns  of  the  table,  the  component  number  and  t-value 
of  a datum.  The  remaining  columns,  following  the  order  of  the  adjust- 
able lambdas,  list  dq,  (t)  / Each  datvim  is  given  in  the  order  in 

which  it  is  stored  internally,  Partials  for  statistical  constraints 
follow  peirtials  for  the  data. 

Option:  column  card  4, 


PARTIAI5  OF  SUM  OF  SQUARES  WITH  RESPECT  TO  ADJUSTABLE  PARAMETERS  - 
TOTAL  SS  - gives  the  partial  of  the  s\im  of  squares  for  all  the 
components  with  respect  to  each  adjustable  parameter. 

FOR  COMPONENT  SS  - gives  the  matrix  of  the  partial  of  the  sum 
of  squares  for  each  component  with  respect  to  each  adjustable  parameter. 
The  rows  correspond  to  the  parameters  and  the  columns  to  the  components. 

Option:  column  5#  card  4, 


PLOT  - a semi -logarithmic  or  arithmeric  plot  of  the  calculated  and 
observed  data.  Each  component  which  has  three  or  more  data  points  is 
plotted  unless  all  points  have  identical  values.  In  the  case  of  the 
semi-log  plot  a)  negative  values  will  result  in  a diagnostic  and  no 
plot,  b)  off  scale  values  are  displaced  by  multiples  of  two  decades  to 
bring  them  onto  the  plot.  If  calculated  and  observed  values  are 
iinequally  cycled,  a C is  printed  at  the  top  of  the  graph.  An  asterisk, 
printed  next  to  the  "t"  value,  indicates  that  the  scale  has  been  stretched 
at  that  point  to  separate  two  data,  with  different  t values,  but  which 
would  have  printed  on  top  of  each  other  because  of  the  limited  resolution. 
An  integer,  printed  at  the  top  of  the  graph,  records  the  ntunber  of  data 
plotted  at  the  same  "t"  value  when  it  exceeds  one. 
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The  symbol  "+'*  in  the  plot  indicates  a calculated  datum.  The 
symbol  in  the  plot  indicates  an  observed  datum.  The  symbol  '*x'* 
indicates  both  calculated  and  observed  data  having  the  same  "print” 
value.  It  should  be  recognized  that  the  plot  has  limited  resolution, 
defined  by  the  separation  of  consecutive  spacings  in  the  printer.  For 
the  semi -log  plot  the  resolution  is  about  5 percent. 

Option:  column  4,  card  4 
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i QO(l)  - represents  q[°(t) 


QC(l)  - represents  q^(t) 


QP(l)  AT  X - is  the  value  that  will  replace  or  be  added  to  the 
current  value  of  F(I,T)  as  instructed  by  the  t-interrupt  control  card, 
when  t = X. 


REORGANIZED  PROBLEM  INFORMATION  - includes: 

1)  the  version  of  SAAM  used,  the  number  of  components  specified 
for  the  model,  and  the  number  of  data  points, 

2)  an  expanded  list  of  the  data  with  modifications  called  for 
by  control  cards  and  with  normalized  statistical  weights, 

5)  the  initial  conditions  of  the  problem,  F(I,0),  changes  in 
F(I)  at  t-interrupts , JUDY(i),  V(I)  and  U(l), 

4)  the  paremeter  values  listed  as  adjustable,  dependent,  fixed. 
Minima  and  maxima  are  given  for  adjustable  parameters.  Relations  of 
dependent  pajrameters  are  expressed  in  terms  of  non-dependent  pao:ameters, 

5)  statistical  constraints  and  f -dei>endence , as  interpreted  by 

SAAM, 

6)  special  output  as  called  for. 


RUNNING  TIME  = X - gives  the  time  (in  computer  clock  units)  it 
took  to  run  the  problem. 

Option : column  1,  caird  1 


SIG  - is  the  mean  wei^ted  variance  of  the  data.  The  initial 
estimate  of  its  value  is  made  from  the  information  in  the  input  deck. 
Subsequent  values  are  based  on  the  solutions. 
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SOUJTION  “ the  first  solution  is  considered  the  zero-th  iteration 
and  is  based  on  the  set  of  parameters  initially  given.  A new  list  of 
psirameter  values  is  given  for  each  t-interrupt  if  the  values  change. 
Included  in  the  output  are: 

1)  the  model  code 

2)  SIG  as  estimated  from  the  input 

5)  X..  values  based  on  the  initial  peirameter  values 
J J 


= -z 


li-)  a table  of  calculated  results  with  headings  for  D ^the  datum 
number),  C (the  component  number),  t,  q.j^(t)//c,  k,  <lj^('t),  Q.j^(t), 

(q^(t)  - q.j5.(t))  and  <lj^(t)/  \{t), 

5)  a weighted  sum  of  squares  for  all  the  data,  for  each  COMP 
and  for  the  block  of  statistical  constraints. 

5)  SIG 


STARTING  TIME  = X - gives  the  clock  reading  at  the  start  of  a 
problem. 

Option:  column  1,  card  1. 


TOTAL  RES  FOR  ITERATION  - the  total  adjustment  to  lambdas  made 
during  am  iteration. 
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PUNCHED  OUTPUT 


A-MATRIX  - the  diagonal  and  upper  half  of  the  A-matrix  for 

adjustable  lambdas  is  punched  out.  The  form  used  is  that 
described  under  matrix  entry  format  in  the  input  section. 
Option:  coliimn  card  4 

COVARIANCE  MATRIX  - the  diagonal  and  upper  half  of  the  covariance 
matrix  for  adjustable  lambdas  are  punched  out  in  matrix 
entry  format. 

Option:  column  2,  card  4. 

CR-VSCTOR  - the  right  side  of  the  noraml  equations  for  adjustable 
lambdas  is  punched  out  as  the  first  row  of  a matrix,  in  the 
matrix  entry  format. 

Option:  column  2 or  5#  card 

MATRIX  CODE  CARD-  one  of  these  cards  is  punched  out  preceding  an 
A-MATRIX  or  COVARIANCE  MATRIX 

column  1 -»■  10  integer  portion  of  problem  number 
column  11  -»•  20  code  for  type  of  matrix 

"1"  for  COVARIANCE  MATRIX 
''3"  for  A-MATRIX 
column  21  -►  35  "the  value  of  SIG 
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i DIAGNOSTICS. 


The  program  contains  a large  number  of  internal  checks  which  may 
resvilt  in  diagnostic  printout.  If  the  problem  calculation  terminates 
because  of  a diagnostic,  the  program  automatically  proceeds  to  the  next 
problem.  This  section  contains  an  alphabetic  listing  of  the  diagnos- 
tics, Those  that  cause  termination  are  marked  with  an  (h).  The  number 
of  the  subroutine  that  produces  the  diagnostic  is  listed  at  the  right, 
Sometimee  a single  input  error  produces  multiple  diagnostics. 
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A SUBSCRIPT  (I)  OF  PARAMETER  (j)  IS  INCORRECT. 


The  subscript  may  not  be  greater  than  the  number 
of  components  in  the  problem. 


ADDING  LAMBDA  (0,j)  HAS  EXCEEDED  LIMIT  OF  I. 

A is  added  to  the  parameter  list  for  evei^ 

component  that  does  not  have  a read  in. 
ADDITIONS  TO  INIT.  COND.  NOT  PERMITTED. 


The  codes  entered  as  initial  conditions  cannot 
be  changed  at  a t-interrupt.  Model  Code  I8. 

ALL  VALUES  IN  COMPARTMENT  I ARE  NEGATIVE.  LOG  PLOT  IS 
IMPOSSIBLE. 

ALL  WEIGHTS  = 1. 

No  weights  were  assigned  to  the  data  and  the 

program  has  assigned  a weight  of  1.  to  each 

datum  and  set  the  ntimber  of  iterations  to  zero, 

ALL  WEIGHTS  FOR  COMPONENT  (l)  ARE  ZERO,  THEREFORE 

KAPPA  (l)  CANNOT  BE  VARIABLE. 

An  adjustable  kapia,  cannot  be  calcxilated  for  a 

component  if  all  data  for  that  component  have 

zero  statistical  weight. 

AMT.  COMP.  I SET  = + 10.E12 

An  amount  calculated  for  component  I has  been  set 
12 

at  + 10  in  order  to  avoid  overflow  or  vmderflow 

on  the  computer  (i.e.  the  calculated  value  was 

12 

outside  the  limit  set  at  + 10  ),  ^todel  codes  1 and  14- 

AT  MATRIX  OVERFLOW. 


The  storage  of  information  has  exceeded  the  available 
storage  locations.  Model  Code 
AT-MATRIX  OVERFLOW. 


Matrices  and  related  information  have  exceeded 
the  available  storage.  Model  Code  11, 
AT-MATRIX  OVERFLOW,  ITER.  SET  TO  ZERO 

Storage  space  for  added  matrix  is  insufficient. 
Only  the  initial  solution  is  made,  and  without 
the  matrix  information. 


(H) 

(H) 

(H) 


(H) 


(H) 


(H) 

(H) 


1 


8 


lii-T 


81^ 

2 


19 


k6 


97 

98 

155 
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AT  T = X EXPONENTIAL  LIMITED  TO  E**50. 

The  amovint  calculated  at  t = X exceeded  e , 
a limit  set  in  the  program.  The  amount  is  set 

50, 

at  e and  the  problem  continued.  Model  Code  9» 

AT  T = X LAMBDA  (l,J)*  Y SET  = I.O69  E 13 

An  exponent,  \t,  has  exceeded  50»  and  been  set 
equal  to  50«  ^todel  code  2. 

AT  T = X LAMBDA  (l,2)  = Y 

The  Gaussian  function  called  for  cannot  be 
calculated  and  f^(t)  for  that  value  of  t 
remains  zero.  Model  Code  I8, 

CODE  INCORRECT  FOR  F(I)o 

The  request  for  calculations  as  entered  under  the 
initial  conditions  heading  incliides  an  incorrect 
comi>onent  number  or  operational  code.  Model  Code  I8. 
COMP  (I)  ISOLATED  FOR  STEADY  STATE  CALC.  SINCE  L(I,I)  = 0. 
COMPONENT  (1)  NOT  IN  MODEL  BUT  INITIAL  CONDITION  GIVEN. 

The  component  named  under  the  initial  conditions 
heading  has  an  index  greater  than  the  hipest 
index  in  the  model, 

COMPONENT  NUMBER  (L)  FOR  DATUM  (l)  EXCEEDS  N. 

A data  point  has  been  entered  for  compartment  L, 

L > N,  the  number  of  components  specified  for 
the  problem, 

COVARIANCE  MATRIX  IRREGULAR 

Negative  diagonal  elements  resulting  from 
inversion  of  ill-conditioned  matrices  may 
cause  this  diagnostic, 

DATA  DEI£TED  FOR  T = X 

The  data  for  t given  by  the  diagnostic  is  incon5>lete 
and  the  point  has  been  dropped  from  the  data  list 
in  the  solution,  ^fc)del  Code  5» 


39 


56 


148 


(H)  147 


(H) 


57 

9 


(H)  2 


41 


68 
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DATA  EXCEED  I. 

I is  the  limit  for  number  of  data  entries, 

DATA  INCOMPLETE 

Some  information  in  the  data  list  is  missing, 

A complete  set  of  observations  must  be 
entered  as  data  at  each  t.  Model  Code  5» 

DATA  INDIV,  I INCOMPLETE, 

Data  giving  the  values  for  individual  I 
are  missing.  Model  Code  11, 

DATA  INSUFFICIENT, 

There  are  not  enough  data  given  (correctly) 
to  continue  the  solution.  Model  Code  5* 

DATA  PLUS  STATISTICAL  CONSTRAINTS  EXCEED  I, 

The  number  of  data  entries  plus  the  number 
of  statistical  constraints  is  limited  to  a 
majcimum  of  I, 

DATA  STORAGE  EXCEEDED, 

Information  is  stored  jointly  vith  the  partials 
matrix.  This  diagnostic  results  when  the  space 
is  insufficient  to  store  all  the  information 
given.  Model  Code 
DEGREES  OF  FREEDOM  LESS  THAN  3, 

The  nximber  of  weighted  data  points  minus  the 
number  of  adjustable  lambdas  is  less  than  5» 

The  solution  will  proceed  if  an  estimate  of  SIG 
exists , 

DEPENDENCE  DIAGNOSTIC.  NO  PARAMETER  (l,J)  IN  LIST, 

The  dependence  relations  include  a parameter 
which  is  not  in  the  parameter  list, 

DEPENDENCE  DIAGNOSTIC,  PARAMETER  (l,j)  IS  NOT  DEPENDENT 
BUT  DEPENDENCE  IS  GIVEN, 

DEPENDENCE  MATRIX  IS  SINGUIAR, 

Insufficient  (or  incorrect)  dependence 
relations  were  entered. 


(H)  2 

68 

(H)  98 

(H)  68 

(H)  7 

(H)  97 

(H)  4 

(H)  6 

(H)  6 

(H)  U 
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DEPENDENCE  ON  PARAMETER  (l,J)  HAS  BEEN  ENTERED  TWICE.  (H)  6 

A single  constraint  has  tvo  terms 
involving  parameter  (l,J). 

DEPENDENT  PARAMETER  (l)  UNDEFINED.  (h)  11 

Give  dependence  relation  for  the  I^  parameter 

DIVIDING  BY  0.  (H)  IO9 

Model  Code  I5. 

ENTRY  CONTROL  CARDS  INCOMPATIBLE.  (h)  2 

Relative  weights  for  some  of  the  data  points 
c€innot  be  combined  with  absolute  weights  for 
other  data  points. 

ERROR  IN  MATRIX  READ-IN.  (h)  155 

The  matrix  supplied  as  additional  information 
Tinder  special  options  (card  2)  has  been 
punched  incorrectly. 

ESTIMATE  OF  ABSOLUTE  DEVIATION  OF  DATA  NEEDED.  (h)  19 

The  degrees  of  freedom  are  less  than  3 and. 
calculation  csuanot  proceed  without 
additional  Information 

EXCESS  DEPENDENTS  =1.  (h)  11 

There  are  I more  dependent  parameters  than 
storage  will  allow, 

F(l)  CODES  ILL-DEFINED.  (h)  14T 

The  codes  entered  under  initial  ccaaditicais 
are  in  question,  ^fodel  Code  18, 

F(l)  = 0.  F(J)  NOT  COMPUTED.  1^4-7 

The  calculation  of  F(j)  is  omitted  (it  is 
set  to  zero)  since  it  requires  division 
by  F(i)  which  is  zero,  ^fc>del  Code  I8, 

FUNCTION  INPUT  ERROR . (h)  97 

Model  Code  3» 

FUNCTION  X NOT  IN  DATA.  (h)  97 

Model  Code  5« 

GIVE  DEPENDENCE  RELATION  FOR  PARAMETER  (l,J).  (h)  6 
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HALT  DIVIDING  BY  ZERO. 

Model  Code  12, 

INITIAL  CONDITIONS  GIVEN  FOR  SUMMING  COMPONENT  (j). 

Initial  conditions  may  not  be  given  for  a summing 
component • 

INITIAL  CONDITIONS  NEEDED. 

See  requirement  under  Model  Code  used. 

INITIAL  ESTIMATE  lAMBDA  (l,J)  OUTSIDE  LIMIT.  ITERATIONS 
SET  TO  ZERO, 

An  upper  or  lover  limit  for  . is  violated  by  the 

^ J 

initial  estimate  of  the  value  of  X. A zero  iter- 

^ J 

ation  calculation  is  performed  and  the  problem  is 
halted. 

JJ  CALCULATION  EXCEEDS  STORAGE. 

The  information  for  X..  calculation  is  stored  in 

J J 

the  dependence  matrix  and  has  exceeded  the  space 
available.  This  can  be  corrected  by  reducing  the 
nvunber  of.  dei>endent  parameters  axid/or,  if  some 
component  numbers  have  been  skipped,  reniombering 
to  reduce  the  total  nvunber  of  components, 

KAPPA  (I ) CANNOT  BE  VARIABLE  BECAUSE  ALL  SIGMAS  INTO 
SUMMER  (1)  ARE  NOT  FIXED. 

Either  k.  or  a.,  must  be  fixed  to  permit  solution 
of  the  model. 

L(I,I)  = 0.  C(l)  ISOLATED  FROM  SYSTEM  FOR  CALC. 

Component  (l)  ceinnot  be  included  for  steady 
state  calculations.  Model  Code  I6, 

LAMBDA  (1,1)  = 0.  COMP.  (l)  ISOLATED  FROM  SYSTEM  FOR 
CALCULATION  OF  U OR  V. 

In  the  calculation  of  steady  state  values 
component  (l)  cannot  be  included  in  the  system, 
LAMBDA  (1,0)  UNDEFINED. 

The  second  subscript  of  a lambda  cannot  be  zero. 


(H)  108 

(H)  5 

(H)  27 

(H)  4 


(H)  11 


(H)  12 


l4l 


42 


(H)  4 
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LAMBDA  MATRIX  SINGUIAR. 

(H) 

141 

The  matrix  of  lambdas  cannot  be  inverted 
Model  Code  I6. 

LAMBDA  (0,1)  NEEDED. 

(H) 

104 

Model  Code  10. 
lAMBDAS  INCCKRECT. 

(H) 

98 

LEunbda  (0,l)  must  be  entered  for  each  component. 
Model  Code  11. 

LIMITS  PARAMETER  (l,J)  INCORRECT. 

(H) 

1 

Lower  limit  higher  than  upper  limit? 
MATRIX  I.D.  FOR  INDIV.  I INCOMPIETE. 

(H) 

98 

The  information  on  the  card  preceding  a matrix 
has  one  or  more  entries  missing.  Model  Code  11. 
MATRIX  INPUT  ERROR. 

(H) 

97 

The  matrix  input  is  incorrectly  punched. 
Model  Code  5» 

MATRIX  INPUT  INCORRECT. 

(H) 

103 

The  matrix  reading  routine  has  found  the  matrix 
dimensions,  as  punched,  incompatible  either 
with  the  dimensions  of  the  matrix  required  by 
the  problem,  or  with  the  subscripts  of  a matrix 
entry,  as  punched. 

MATRIX  N.G.  FOR  I. 

(H) 

98 

The  matrix  input  for  individual  I,  is 
incorrectly  punched.  Model  Code  11. 
MATRIX  SINGUIAR. 

(H) 

57 

The  matrix  used  in  the  solution  is  singular. 
Model  Code  3» 

MCDEL  NOT  SUITABLE  FOR  ANALYTIC  SOUJTION. 

(H) 

17 

The  model  requires  the  use  of  the 
differential  equations  solution,  model  types  1 or 
NO  ADJUSTABLE  lAMBDAS.  ITERATIONS  SET  TO  ZERO. 

If  there  are  no  adjustable  lambdas  the  number 
of  iterations  is  set  to  zero. 

NO  ADJUSTABIE  PARAMETERS. 


67 
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NO  PARAMETER  (l,j)  IN  LIST  BUT  STATISTICAL  CONSTRAINT  GIVEN. 
Ktrameter  (l,J^  entered  in  a statistical 
constraint, is  not  in  the  list  of  i«xameters, 

NO.  OF  ADJUSTABLE  PARAMETERS  EXCEEDS  I. 

NO.  OF  DEPENDENTS  EXCEEDS  I. 

NO.  OF  PARAMETERS  EXCEEDS  I. 

NO.  OF  STEPS  NEG.  FOR  CALC.  OF  DATUM  I. 

An  inconsistent  t-interrupt  entry  has  probably 
resulted  in  a negative  t increment. 

OVERFLOW  COMPUTING  F(I)  AT  T = X. 

The  solution  in  progress  is  completed  and  the 
problem  is  terminated,  ffodel  Code  18. 

OVERFLOW  IN  COMP.  I AT  T = X. 

Model  Code  I5. 

PARAMETER  DEPENDENCE  ON  C(l)  INCORRECT. 

I probably  exceeds  the  highest  component  index. 
PARAMETERS  I AND  J ARE  DUPLICATES. 

A peirameter  has  been  entered  twice  in  the 
input  data. 

PARAMETER  I = X,  NOT  PUNCHED. 

The  limits  of  the  floating  point  field  were 
insufficient  to  accommodate  the  parameter  value. 
PARTITION  I OF  THE  MATRIX  IS  SINGULAR.  THE  FOLLOWING 
ADJUSTABLE  PARAMETERS  ARE  INVOLVED. 

Symmetric  matrices  are  peurtitioned  into  independent 
blocks  for  inversion.  A partition  involving  the 
listed  parameters  is  singular.  At  least  one  of  the 
parameters  in  this  block  is  dependent. 

RATIO  OF  CALCULATED  SIGMA  TO  ESTIMATED  SIG  = X,  ITERATIONS 
SET  TO  ZERO 

The  deviation  of  the  calculated  from  the  observed 
values  is  more  than  ten  times  greater  than  implied 
by  the  input  estimate  of  the  errors.  Better  initial 
estimates  of  iwirameters  or  larger  estimates  of 
standard  deviations  of  data  are  needed. 


(H) 


(H) 

(H) 

(H) 

(H) 


(H) 

(H) 

(H) 


(H) 


5 

1 

1 

21 


114.7 


109 


19 


156 


29 
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SIG  = 0.  RES.  NOT  MODIFIED 

The  RES  vector  resulting  from  solution  of  the 
normal  equations  could  not  be  modified  because 
SIG  =0,,  and  modification  includes  a calcula- 
tion requiring  division  by  SIG.  The  solution 
continues  with  the  unmodified  RES. 

SIGMA  - KAPPA  MATRIX  IS  SINGULAR. 

The  matrix  of  coefficients  of  the  normal  equations 
used  to  calculate  <7*s  and  k*s  is  singular 
SIGMA  KAPPA  MATRIX  SINGULAR. 

The  matrix  of  normal  equations  for  the  solution 
of  sigmas  and/or  kappas  is  singular. 

SPREAD  OF  VALUES  FOR  COMPARTMENT  I EQLJAIS  ZERO. 

Plot  will  be  omitted. 

SQUARE  ROOT  UNDEFINED. 

The  result  of  the  calculation  would  be  an 
imaginary  number.  Model  Code  12. 

STD.  DEVS.  FOR  I GIVEN  WITH  MATRIX  FOR  J. 

STORAGE  OVERFLOW. 

The  available  storage  is  insufficient  for  the 
calculation  of  codes  submitted.  It  may  be 
possible  to  simplify  the  input  to  fit  within 
available  storage.  Model  Code  18, 


25 


(H)  54 


84 

(H)  108 


(H)  98 

(H)  147 


T,  AT  T-INTERRUPT  I EXCEEDED  BY  DATUM  J.  (H)  2 

The  T value  of  datum  (j),  in  the  block  terminated 
by  T-inte’brupt  card  (l),  exceeds  the  T value  of 
the  T-interrupt. 

T-INTERRUPT  DIAGNOSTIC.  A CCMPONENT  NO.  EXCEEDS  N.  (H)  52 

The  subscript  of  a t-interrupt  entry  exceeds  the 
number  of  components  in  the  model. 

T-INTERRUPT  DIAGNOSTIC,  PARAMETER  (l,J)  OR  PARAMETER  (K,L) 


NOT  IN  LIST. 

A parameter  to  be  changed  at  a t-interrupt  is  not 
in  the  original  parameter  list. 


(H)  52 
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T-INTERRUPTS  EXCEED  I.  (H)  2 

I is  the  maximum  nvunber  of  t-interrupts  permitted, 

THE  LAMBDA  MATRIX  IS  SINGULAR  4-2 

In  the  calculation  of  compartment  sizes,  V,  the 
X matrix  is  singular.  Model  Code  1, 

TIME  ALLOWED  EXCEEDED  AT  T = XX,  INCREASE  TIME  FACTOR 

IF  DESIRED,  (h)  21 

Program  estimates  computation  time  and  limits 
computations  to  a preset  internal  reference.  This 
can  be  changed  by  a Time  Factor  entry  on  card  5» 

TOO  FEW  STEADY  STATE  EQUATIONS,  (H)  97 

The  number  of  equations  must  be  equal  to  the 
n\amber  of  components.  Model  Code  3» 

TOO  MANY  COMPONENTS  IN  MODEL.  (H)  2 

TOO  MANY  KAPPAS  LISTED.  (h)  3 

TOO  MANY  STEADY  STATE  EQUATIONS.  (h)  97 

The  number  of  equations  must  be  equal  to  the 
number  of  components.  Model  Code  3« 

V AND  U MAY  NOT  BOTH  BE  GIVEN.  42 

No  compartment  sizes  (v)  may  be  ^ven  if  any 
steady  state  inflow  rates  (u)  are  given,  and 
vice  versa.  Model  Code  1, 

VALID  ONLY  FOR  NON-SUMMERS.  l4 

Information  already  available  in  storage  is 
printed  out  but  is  valid  only  for  the  f^  of 
the  model, 

VARIANCE  NEGATIVE  69 

A stsindard  deviation  cannot  be  calculated  because  a 
diagonal  element  of  the  covariance  matrix  is 
negative . 

VARIAJICE  NEGATIVE,  S,  D.  SET  TO  ZERO  42 

A diagonal  element  of  the  covariance  matrix  is 
negative  and  its  square  root  (the  standard 
deviation)  cannot  be  calculated. 
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WEIGHTS  FOR  COMPONENT  (l)  ARE  ZERO  THEREFORE 

SIGMA  (J,I)  CANNOT  BE  ADJUSTABIE.  (h)  19 

An  adjustable  sigma  cannot  be  calculated  If  the 
data  for  the  svunmer  component  have  zero  statistical 
weight. 

ZERO  DIAGONAL  ELEMENTS  MATRIX  IS  SINGULAR  (h)  kk 

One  or  more  diagonal  elements  of  the  symmetric 
matrix  being  inverted  are  zero, 

ZERO  DIVISION  IN  COMP.  I.  T = X 45 

Calculation  of  the  I^^  component  called  for 
division  by  zero  at  the  point  whose  T value 
is  X.  Model  Code  7* 
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j Model  Types 

I 1.  Linear  differential  equatloiis  with  constant  coefficients; 


'if.('t’)  A ip 

^ = L - L 

i=l  i=0 

i/j  i/j 

Initial  conditions : f . (O) 

J 

jA 


The  are  time  independent.  Y.(t)  is  an  arbitrary  input  function 

ij  J 

into  Equivalent  schematic 


Computational  procedure, 


th 


The  differential  equations  are  solved  using  a 4 — order 


Runge-Kutta  method 
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In  addition  and  independently  of  the  above  solution,  a steady 
state  solution  V can  also  be  obtained  from  a constant  injut  vector  U, 
entered  separately.  In  aiatrix  notation  this  is  given  as 

V = u 

>.  / ’ 

or 

u = X V 

where  U and  V are  vectors  and  X is  the  matrix  of  coefficients : 


^11 

-X12 

■^15 

-X21 

X22 

-^23 

-Si 

"X32 

S3 

• 

# 

• 

« 

• 

At  present  these  calculated  values  cannot  be  fitted  against 
observed  values. 

Special  Inputs : 

Card  4:  enter  "1"  under  MODEL  CODE.  This  entry  may  be  left 

blank,  in  which  case  the  program  will  try  to  use  model 
code  2,  if  the  problem  meets  the  requirements  for  that 
model  code.  If  not,  it  will  use  model  code  1, 

Initial  conditions:  enter  initial  conditions  for  differential 

equations,  f^(0),  under  INITIAL  CONDITIONS. 

Enter  either  the  vector  V or  the  vector  U in  their 
appropriate  fields,  if  a steady  state  solution  is 
desired.  VHien  an  entry  (or  V^)  is  omitted  it  is 
autcmatically  set  to  zero  by  the  program. 


VI  - 5 


August  1966 


Special  Outputs : 

INVERTED  lAMBDA  MATRIX  (Option). 

U(l),  V(l)  smd  R(I,J),  when  steady  state  solutions  are 
requested. 

R(I,J)  s X(I,J)*V(J) 
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2.  Sums  of  Exponentials 


= * 


«k  = 'Sl  X '’■=1  '*1  ' 


Equivalent  schematic 


Computational  procedure : 

~\o.t 

Evaluation  of  e for  each  datum 

Model  code:  2 (column  50>  card  4) 

User  may  leave  code  "blank"  in  which  case  program  will  choose 
model  code  2 if  it  can,  or  else  it  will  choose  either  model  code 
1 or  4,  whichever  is  appropriate. 
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Similarity  Transformation  (MAPPING) 

This  SOLVE  routine  i)erforms  a similarity  transformation  on  a given 
matrix  M: 

F = XM\"^ 

F,  X and  M are  rxn  matrices. 

Option  A: 

1 = 0,  1,  . . n 
J — 1,  2,  • • . , n 


fl(t.j)  = F(I,J) 


where  F(l,j)  is  the  I row  and  J column  element  of  the 


matrix  F for  I / 0,  and  F(0,j)  = ^ F(I, 

1=^1 


7) 


= U(I) 

= V(I) 

where  U(l)  and  V(l)  are  solutions  of  the  simultaneous  set  of 

2n  linear  equations 
n 


^ F(I,J)V(J)  - 11(1)  = 0 


X — n. 


= C(I) 


y A(I,K)V(K)  + z B(I,M)  U(m)  ^ I = 1,  2,  n 

ifel  Mel 

A(I,K),  B(I,M)  and  C(l)  are  required  for  the  calculation  of 
U(l)  and  V(l)  and  are  entered  by  the  user  under  special  input 


described  later 
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= R(I,J) 


vhere  r(I,J)  = fi(t  )xv(j).  The  r(I,J)  are  calculated 
^ J 

after  V(j)  are  determined, 

= X(I) 

where  X(l)  = ^ D(I,K)v(k)  + ^ E(I,J)u(j) 

K J 


D(I,K)  and  E(I,J)  are  entered  as  special  input  described 
later.  The  X(l)  are  calculated  after  the  V(K)  and  U(j)  eu:e 
determined. 


Option  B 


= - F(I,J) 

for  I J and 

0 

M 

= F(I,J) 

for  I = J and 

1 = 0 

f same  as  optical  A 

J 

= fi(t^j)v(j) 

same  as  option  A 

(Note:  In  linear  compartmental  systems  with  constant  coefficients  U, 

V,  R and  X are  steady  state  quantities). 

Special  inputs : 

Card  2:  "n"  (dimension  of  F(I,J))  under  NO  OF  COMPONENTS, 

Card  4:  ”3"  under  MODEL  CODE 


Column  70 


f ”1"  Option  A 

V Blank  Option  B 
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! Data: 


COMP  NO. : "1."  for  fi 


"2.'*  for  fa 


T:  I + .01*J  for  F(I,J)  [e.g.  1.02  for  F(1,2)  1 

50.  + .01*1  for  V(l)  [e.g.  50.11  for  V(ll)  3 

60.  + .01*1  for  U(l)  [e.g  60.02  for  U(2)  1 

100.  + I + .01*J  for  F(I,J)  [e.g.  101.02  for  R(1,2)  1 

150.  + .01*1  for  X(I)  [e.g.  150.12  X(12)  1 


Initial  Conditions : 

No  entries  needed  except  when  fa(t)  computations  are  called  for, 
in  which  case  enter  a "1."  in  columns  {h-2  -►55)  of  the  '*26** 
termination  card  for  this  block. 

Lambdas ; 

Elements  of  \(l,j)  matrix 
Special : 

l^trix  M is  entered  after  last  "26"  termination  card  of  regular 
data  deck.  The  format  is  as  given  in  Section  III  28.  Each 
matrix  is  terminated  by  a card  with  "26"  in  column  5 
The  constants  A(I,K),  B(I,M)  and  C(l)  involved  in  the  equations  for 
U and  V: 


+ 


C(I) 


K 


M 


are  entered  after  the  matrix  M block,  one  card  for  each  term  of 
the  equation  in  the  following  format: 
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column  4 -►  5 enter  30  for  V(k),  60  for  u(m)  (integer) 

column  9 10  enter  K or  M (integer) 

colximn  13  -*•  25  enter  A(I,K)  or  B(I,M)  (decimal) 

On  the  card  containing  the  last  term  of  each  equation  make  the 
additional  entries: 

column  26  '*=” 

column  27  40  enter  C(l)  (decimal) 

The  program  requires  that  n equations  be  entered.  The  entire  set 
of  equations  is  terminated  by  a card  with  a "26"  in  column  4 ■>  5<> 
The  constants  D(I,K)  and  E(I,J)  involved  in  the  equations  for  X: 


X(I) 


E(I,J)  U(J) 


are  entered  next  under  the  following  format,  one  card  for  each 
term: 


column  5 -*•  5 enter  "15O 
column  9 -*■  10  enter  I 
colvmm  19  ■*  20  "30"  for  V(k), 

column  24  -*■  25  enter  K or  J 
colxomn  27  40  enter  D(l,j)  or 


(integer) 

(integer) 

"60"  for  U(j)  (integer) 
(integer) 

E(I,J)  (decimal) 


The  entire  set  of  X(l)  entries  are  terminated  by  a "26"  in  column 
4 -*•  5«  The  x(l)  values  are  entered  under  Data  as  described  earlier 
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j 

i 

Special  outputs 

! 

j 

J lAMBDA  MAlTRIX  - is  the  matrix  X to  be  used  to  solve  the  equations 

j FV  = U. 

j LAMBDA- INVERSE  MATRIX  - is  the  inverse  of  the  lambda  matrix. 

MATRIX  - is  the  matrix,  M,  given  in  the  input  for  solution  of  the 
equations  F = X M X”^. 

STEADY  STATE  EQUATIONS  - are  linear  equations  in  U and  V to  be 
solved  given  the  matrix  of  X*s 

STEADY  STATE  FUNCTIONS  - are  additional  constraints  in  U and  V. 
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4.  Ordinary  non-linear  differential  equations 
n n 


df . 
J 


dt 


= Z Si  - Z Sj 


i=l 


1=0 


j = 


initial  conditions : f . (O) 

J 


%Jt)  = \ ).  fj(t)  or  q.(t)  = K-  f.(t) 


q (t ) = K ‘ ) cr.  . f . ("t ) 


The  parameters  designated  hy  a 'prime*  can  he  f-dependent. 
For  exaunple,  let 
df. 


dt 


- “ ^21^1  ^51^1  ^15^3 


If  ^21  “ ^21/^ j-’  ^31  ” ^31^^i^  ^13  ^13^ 


^13  '^13^ 


df. 


dt 


= i-iaVa  - ^21  Vs  ^ hi  ^ VS 


Similarly,  if  q^  = ^7  “ ^7/^6^  then  q^  = K^f^/f^ 


Y.(t)  is  an  arbitrary  input  function  into  j. 


' } Up 


Coraputational  procedure : 

4th  order  Runge-Kutta  method  is  used  for  solution  of  differential 

I 

j equations . 
j Special  inputs : 

Card  4 - enter  "4"  under  MODEL  CODE.  If  this  entry  is  blank 

the  program  ■\rill  appropriately  choose  bet^reen  models  2, 
or  4. 

\’s  ct's,  k’s — enter  f-dependence  information  in  columns  57  59* 

Special  outputs ; 

In  case  of  interruptions  in  solution  special  intermediate  results 


are  printed 
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5.  Linear  Combination  of  Spectra 


3 k 


r- 


j - 2, 5>  • • • 


where  f.(t,  ) are  given  values  for  spectrum  j at  t,  and 

J j£  K 

is  some  arbitrary  t shift  in  the  coordinate  system  for  f.. 


II 


jA 


Computational  procedure:  linear  regression 

Special  inputs : 

Card  k:  enter  "5'*  under  MODEL  CODE 


Data 


COMPONENT  - index  j of  f. 


k = 2,5,. 


T - value  of  t,  . Every  function  f.(t,  + X.,)  must 
■K  J J ^ 


have  a value  at  each  t,  . 

k 


OBSERVED 


VALUE  - f.(t,  + X..)  and  q.v('tv)  values 

,]  K -IX  K K 


9 - derivative  of  f . (t,  + X . , 

J k 


if  available 
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6 . Sum  of  Gauss ians  and  Sxponentials 


f,(t)  = 

J 


4^  X 


js 


= e 


i V 

4l^(t)  = Kj.  ) 

ilk 


for  3 > 3 


for  h > 5 


i > 5 
k > 3 


Special  input : 

Card  14-  - "6"  under  MODEL  TYPE 


I 
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7.  Population  Survival  P^Jnction 


“X.-|t  -X.pt  X., 

f.(t)  = f.(0)  e [1  - (1  - e 

J J 


Special  Input  ; 

Card  4 - "7"  under  MODEL  CODE 


for  j > 3 


for  k > 3 
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8.  Power  Series 


f.(t)  = 

^ for  j = 1,  2,  ...,m 

= 1. 

m+1 

= 

jA 

Special  inputs 
Card  4 - 

"8"  under  MODEL  CODE 

"m"  in  column  70 
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9.  Special  Function 


fj(t,s)  = (1  + Xj^il)e 


= \ > ''kj 


Special  Inputs : 

Card  k - "9"  under  MODEL  CODE 


j>  3 
k > 3 
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j 11 . Population  Mean  and  Covariance  Matrix 

j Purpose : Given  a population  of  k studies  in  which  each  study  ^ 

! contains  a set  of  n parameters  with  values  X.  (X.^,  X.„,  ...  X.  ) and 

I — 1 ' il^  i2^  in' 

j covariance  matrix  find  the  mean  pareimeter  values  X (X^,  Xg,  ...  X^) 
and  the  covariance  matrix  (v)  for  the  population. 

Solution  equations  ; the  program  solves  the  following  set  of 
simultaneous  equations,  approximate ly 


"1 

"1 

- * 

X = 

^ (V.  + V) 

) (V.  + V)  X. 

1 

Im 


1 

lc-1 


. „ (x.  - Xj  (X.  - X ) 

1^  im  ' li  i'  ^ im  m' 


where 


element  of  covariance  matrix  V 
4 m 


statistical  "jeight  of  parameter  study  1,  internally 


calculated. 

Special  inputs: 

Card  2 - irGMSSR  OF  COlTOh’ENTS  = n = number  of  parameters 


IJUIC3ER  OF  ITSRAliONS  ; Blank 

Card  -I- 

ilODEL  CODE;  "ll" 

/-I  \ 

Data'^“'^  : 

COMP:  Parameter  mimber  (say,  j_) 

f;  studj'-  number ^ (any  number  0.  ->  999.) 

OBSERVED  VALUE : the  value  of  X. . - parameter  ^ in  study  i 

J 
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IISIGHP : enter  "1  directly  or  through  control  card 

Lambdas ; 

COMP  IIUMBER;  subscripts  for  entered  for  j = 1^  o..^  n. 
IJo  other  entries  are  required  under  this  format. 

Additional  Input : 

Immediate!;'-  after  the  standard  deck  the  follo-^/ing  are  entered  in 
the  order  indicated. 

For  each  study ' 

Control  card; 

Columns  1 ->  10 ; The  study  nunfoer  ^ (integer) 

11  ->  20:  Code  (integer)  for  the  type  of  matrix  input 

which  is  to  follow 
"l"  - covariance  matrix 

”2"  - correlation  coefficient  matrix  with 
standard  deviations  on  diagonal. 

"3"  - normal  equations  matrix 
"i]-"  - normal  equations  raatrix  follov^ed  by 
vector  of  standard  deviations 
(Any  one  of  the  above  may  be  used  as  input 
matrix) 

21  ->  55 : SIG  for  study  ^ (E  format ) . 

(SIG  = s-um  of  squares  of  deviations  dl-vlded 
by  degrees  of  freedom.) 

]vfe,trix^^^ : See  Section  111-28 


VI  18 


Aug,  1966 


Vector:  'Seune  as  matrix  format  (Section  III- 2.1) 

'This  entry  is  required  only  if  called  for  by  code 
entered  on  "control  card". 


Footnotes : 

(1)  Punched  cards  for  these  entries  can  be  obtained  in  proper  format 
for  each  study  i by  entering  a "l"  -under  options  in  either 
columns  5 or  5^  when  "running"  study  ^ under  its  model  code. 

(2)  Study  numbers  must  be  different  from,  each  other  in  the  integer 
portion-i  They  are  ordered  internally  by  their  integer  values 
and  assigned  internal  n-umbers  from  1 ->  k. 

(5)  Since  the  matrices  are  symmetric,  only  upper-right  or  lovrer  left 
portions  (including  diagonal)  need  be  entered. 
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Special  outputs 


A-MATRIX  FOR  POPUIATION 

A-MATRIX  (optional)  is  the  A-tfetrix  for  an  individual  in  the 
population. 

Option:  Column  6,  Card  ^4-. 

CORRELATION  COEFFICIENT  MATRIX  FOR  POPULATION 

INDIVIDUAL  =1,  CODE  = J,  SIG  = X - is  the  code  card  identifying, 
for  the  matrix  which  follows,  the  problem  (INDIVIDUAL),  the 
type  of  matrix  (CCOE)  and  SIG.  The  Code  is  "l'*  for  a covari- 
ance matrix,  "2”  for  a correlation  coefficient  matrix  with 
the  standard  deviations  as  dieigonal  elements,  ”3”  for  the 
A-matrix  and  ”4"  for  the  A-matrix  followed  by  the  vector 
of  standard  deviations. 

MATRIX  AS  READ  - is  the  matid.x  supplied  to  the  program  for  an 
individual  in  the  population  being  studied, 

SIG  FOR  POPULATION  = X - is  the  mean  variance  of  the  parameters 
for  the  population  being  studied. 

VARIANCE  CO-VARIANCE  MATRIX  FOR  ADJ.  POPULATION  - is  the  covari- 
ance matrix  for  the  population  calculated  after  the  param- 
eters (for  each  individual  in  the  population)  have  been 
adjusted  toward  the  mean  within  a 95^  confidence  limit. 


SECTION  VII 


COMRJTATIONAL  PROCEDURES  AND  METHODS 
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The  followxng  are  the  computational  stages  in  SAAM  connected  with  the 
solution  of  a problem: 

Read-in 

Set-up 

Solve  (zeroeth  iteration) 

Iterate:  a)  ftirtial  derivatives 

b ) Parameter  adjustment 
Final  Results 
Wind-up 
Special 
Plot 

READ-IN 

A problem  deck  consists  of  a nximber  of  distinct  blocks,  each  carry- 
ing a different  type  of  information  about  a problem.  Some  computations  are 
intermixed  with  read -write  instructions,  but  most  occur  after  the  deck  is 
i read  in. 

The  first  four  caj:ds  of  the  problem  deck  are  stored  in  core.  When 
"clocking”  is  requested,  the  clock  time  at  the  start  of  a problem  is 
stored  sind  printed. 

The  DATA  block  is  read  and  stored  in  core  next.  Instructions  carried 
by  the  Data  Control  Cards  eire  executed  diiring  the  read -in. 

INITIAL  CONDITIONS  and  the  KAPPAS,  LAMBDAS,  and  SIGMAS  blocks  are 
read  after  the  DATA. 

DEPENDENCE  RELATIONS  are  read  into  temporary  storage.  Dependent 
parameters  are  expressed  in  temns  of  independent  ones  (adjustable  and 
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fixed) : 


Ax  = By 


where  x is  the  vector  of  the  dependent  parameters  and  y is  the  vector  of 
the  independent  ones.  A and  B are  matrices  of  coefficients.  Solution 
for  X gives: 


X * a“^  By 

The  matrix  A B is  stored  for  suhseq.uent  use. 

STATISTICAL  CONSTRAINTS  are  read,  and  stored  jointly  with  the  data 
to  be  used  in  the  least  squares  fit. 

T-INTERRUPT  CHANGES  IN  f^ (t ) are  read  and  stored  in  core . 

T-INTERRUPT  CHANGES  IN  PARAMETERS  are  stored  on  tape  and  brought 
into  core  each  time  they  are  needed.  | 

Inputs  required  for  special  model  types  are  read  during  the  set-up 
wind-up  and  special  computational  stages  of  the  problem  solution.  ^ 


Reorganization  of  data  and  peirameters  takes  place  during  and  after  ' 
the  read  in.  Data  are  rearranged,  new  parameters  are  added  when  required,  ? 
internal  codes  are  set  up  for  computational  control,  etc.  Some  of  the  j 
reorganized  information  apj)ears  in  the  print-out. 


The  input  DATA  are  rearranged  in  order  of  increasing  t within  each  ' 

_ , 

t-interrupt  block.  In  some  model  types  data  are  added  or  deleted.  A \ 

I 

diagnostic  usually  explains  deletions.  j 

Statistical  weights  are  assigned  to  the  data.  First,  weights  are  j 

'J 

assigned  in  accordance  with  the  weight  code.  When  wei^ts  are  assigned  to  5 

I 

the  data  directly,  they  are  entered  as  is.  When  a standard  deviation  (s.d.)> 
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I 

; is  assigned  to  a datum  its  weight  (W)  is  calculated  as: 

i 

I 

I W = 

I (s.d. 

j 

I When  a fractional  deviation  (f.d.)  is  assigned  to  a datum,  the  weight  (w) 

i 

is  calculated  as 

j W = i 

j [(f.d.)  * q°(t)1® 

I 

< O / \ 

j where  q (t)  is  the  OBSERVED  value  for  the  datiim.  When  the  standard 
deviation  (s.d.)  or  observed  value  q°(t)  are  zero,  the  weight  assigned 
to  the  datum  is  also  zero. 

After  weights  are  assigned  to  all  data,  they  are  normalized  so  that 
the  s\am  of  the  weights  equals  the  number  of  data  points  having  non-zero 
weights . 

Weights  assigned  directly  to  data  are  considered  relative,  whereas 
weights  derived  frcM  standard  or  fractional  deviations  are  considered 
absolute . The  program  does  not  permit  the  mixing  of  relative  and  abso- 
lute weight  assignments  in  the  input  data. 

When  the  wei^ts  are  absolute,  it  is  possible  to  calculate  a '‘read- 
in”  average  variance  (SIG)  per  datum: 

SIG  , . = (W  )-  * (s.d.)  2 

read-in  ' n'k  ‘ 

where  (W  ),  is  the  normalized  wei^t  for  any  datum,  k,  having  non-zero 
H i£ 

weight, 

Wei^ts  (W^)  are  assigned  to  the  statistical  constraints: 

SIG  , . 

,,  read-in 

^ " (s.d.)/ 


Vll-k 


August  1966 


where  (s.d.)  is  the  standard  deviation  read  in  for  the  statistical  con- 
straint , 

SET-UP 

Computations  to  test  and  augment  special  information  required  for 
the  particular  model  type  are  carried  out  here.  Every  model  type  has 
its  own  set-up  routine. 

SOLVE  (Zeroeth  Iteration) 

Ctoce  the  parameters,  data  and  other  entries  are  organized,  a 
solution  for  the  particular  model  type  in  the  problem  starts. 

When  a solution  starts  with  a t-interrupt  all  changes  in  f^  and 
parameters  are  executed  first . Values  for  the  dependent  parameters  are 
then  calculated  from  the  stored  dependence  relations. 

The  solution  proceeds  in  accordance  with  the  equations  defining  the 
model  type,  fj^(t)  values  are  derived  for  each  datum  and  for  each  of  the 
statistical  constraints. 

Sigmas  (and  kappas)  to  yield  a least  squares  fit  of  the  data  are 
calculated  next  from  the  linear  regression  equations : 


for  all  k. 


The  calculated  cr^^,  are  then  used  to  calculate  qj^(t).  Comparison 
of  and  permits  the  calculation  of  a 


SIG 


calc 


y ^(t)  - i^(t)r  \ 

k 


d.f. 


where  d.f.  are  the  degrees  of  freedom,  the  nximber  of  weighted  data 
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points  less  the  number  of  adjustable  lambdas. 


In  subsequent  computations,  is  used  in  place  of  SIG^ 


read-in 


if  it  value  is  lower, 

Wei^ts  for  statistical  constraints  are  recalculated  when  the  value 

of  SIG  - is  lower  than  the  value  of  SIG  , . : 
calc  read-in 


after  statistical  wei^ts  are  reassigned,  wei^ted  sum  of  squares  and 
SIG  values  are  calculated  for  all  data  and  statistical  constraints. 


100,  the  problem  solution  is  terminated  at  the  zeroeth  iteration,  imply- 
ing that  initial  estimates  for  parameters  are  not  good  enough  to  continue. 
This  teiTTii nation  cem  be  overcome  by  improved  initial  estimates  of  param- 
eters or  by  entering  larger  steuadard  deviations  for  data,  (This  will 
not  alter  final  least  squares  parameter  values  or  their  estimated 
uncertainties  if  statistical  constraints  are  altered  by  the  same  factor)*. 
The  latter  change  should  only  be  made  if  convergence  with  poor  initial 
estimates  is  jxDssible  or  economical. 


W 


c 


SIG^g^^^  is  compared  with  SIG 


read-in  * 


If  the  ratio  is  greater  than 


*Tentatively  this  can  also  be  overcaae  by  entering  a large  value  of  P 
on  card  5,  P = 5.  is  routinely  used  internally.  P = 10000,  would  in 
almost  every  case  evade  the  test. 
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ITERATIONS 


A.  I^rtial  Derivatives 


After  the  zeroeth  iteration,  partial  derivatives  of  the  calcu- 


lated <Ij^(t)  with  respect  to  each  of  the  adjustable  lambdas  are  derived. 


The  calculations  are  performed  by  numerical  approximation 


^\(t) 


Aqjj.(t) 


ax. . 

ij 


^^ij 


Each  of  the  adjustable  X.,  is  changed  by  a small  fraction,  P^,  (Card  5) 


to  a value  (X.  . + Pi  * X.  .)  and  a new  q,  (t),  designated  as  q?(t)  is 
ij  ij  K.  K 


obtained  using  the  procedure  outlined  under  SOLVE.  An  approximation  for 
the  partials  is  obtained  as 


- q^(t) 


* ^ij 


The  value  of  P^,  normally  .01,  can  be  changed  on  Card  The  recalculated 
q^(t)  include  readjusted  sigmas  and  kappas,  when  the  latter  are  adjustable 
or  dependent.  lartial  derivatives  are  calculated  for  all  data  and 
statistical  constraints. 

Normal  equations  are  derived  from  the  matrix  of  partial  derivatives, 
and  an  estimate  for  a correction  vector  (res)  for  the  adjustable  X^ 
based  on  linear  approximation  theory  — is  calculated: 


RES  = (a'^Wa)"^  a^Wb 


a is  the  matrix  of  partial  derivatives,  W is  a diagonal  matrix  of  the 


weights  and  b is  the  vector  [qj^(’t)  “ qj^Ct)"!. 
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! B.  Parameter  Adjustment 

I After  the  correction  vector  is  calculated  the  RES(l)  are 

j further  modified  one  at  a time,  by  a factor 

' 2 

' [P  * RES(I)  ] 

j VAR(I) 

I 1 + [P  * RES(I)  f 

j VAR(l) 

I vhere  P is  a reeid-in  constant  (Card  5)  and  VAR(i)  is  the  calculated 

I 

i variance  for  variable  I,  (When  P is  not  entered  a value  of  5»  is  auto- 

I 

matically  assigned  to  it.)  Basically,  when  a variance  is  very  large 
compared  to  the  calculated  (RES)^,  the  factor  is  very  small,  and  when 
the  variance  is  small,  the  factor  goes  to  unity.  Modification  of  the 
RES(l)  starts  with  the  last  adjustable  one.  Once  it  is  adjusted  the 
remaining  unmodified  RES(l)  aure  recalculated  subject  to  acceptance  of 
the  values  for  the  modified  ones.  Variances  for  the  unmodified  RES(l) 
are  recalculated  and  used  for  subsequent  adjustments. 

The  modified  RES  vector  is  added  to  the  vector  of  adjustable  X*s, 
and  the  values  are  tested  against  the  inq>osed  upper  and  lower  limits. 

If  a limit  is  violated  the  total  correction  vector  RES  is  scaled  down 
to  stay  within  the  most  restricting  limit.  The  scaling  factor  is 
called  CONAB. 

A new  solution,  Q.^('t)  is  obtained  for  the  adjusted  lambdas.  These, 

together  with  the  %.{%)  of  the  previous  solution  and  the  observed  values 
£ 

q°(t)  are  then  used  to  scale  the  RES  correct!  C3n  vector  farther  to  obtain 

lx 

an  improvement  in  the  fit  of  the  data.  This  scaling  is  based  on  a linear 
extrapolation  using  the  relation: 
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%i^)  - ‘luC't) 

CONj^  = — — k = 1,  m 

- q^(t) 

from  which  an  average  scaling  factor,  CON,  is  calculated.  Because  of 
non-linearities,  this  adjustment  is  repeated  until  the  new  sum  of 
squares  improves  by  less  than  (l-CONMIN)  of  its  previous  value — or  after 
three  tries  in  any  one  iteration.  The  value  of  CONMIN  is  read  in  on 
Cand 

If  the  svim  of  squares  becomes  worse,  the  RES  vector  is  cut  down  by 
a factor  of  ten,  and  the  procedure  is  repeated. 

The  best  sum  of  squares  obtained  during  the  various  adjustment 
stages  is  recovered  as  the  final  solution  of  the  xteration.  A single 
iteration  involves  the  calculation  of  partials  and  the  subsequent  sd- 
justment  of  the  correction  vector. 

FINAL  RESULTS 

A problem  terminates  its  iterations  when  the  maximum  number 
specified  (Card  2)  is  reached  or  when  fractional  improvements  in  the 
sum  of  squares  is  less  than  (1-E).  The  value  of  E is  entered  on  Card 
The  best  solution  is  retrieved  for  the  final  results. 


SECTION  VIII 


PROGRAM  ORGANIZATION 


V 
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The  control  of  the  various  canputational  stages  of  SAAM  is 
governed  hy  a single  main  program  in  an  overlay  version  (FORTRAN  IV) 
amd  by  a series  of  main  programs  in  a chained  version  (FORTRAN  II), 
Schematically  the  two  versions  may  be  described  as  follows : 

OVERLAY 


CHAIN 


Maltiple  links,  each  governed  by  MP5  and  containing  selected 
solution  routines  may  be  set  up  for  the  CHAIN.  Dummy  routines  are  used 
as  needed  to  minimize  core  req\iirements. 

MAIN  roOGRAM  (for  overlay) 

A single  main  program  governs  the  program  flow. 
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MIN  PROGRAMS  (for  CHAINS) 
MAIN  100 

MP  0 


MIN  1 


MIN  2 


MAIN  5 


MAIN  30 


MIN 


MP  1 


MP  2 


MP  5 


MP  50 


1^5 

51 

52 
58 
80 


Iff  51 


89 
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READ-IN 


READO 

1 

READl 

2 

READ2 

3 

READ5 

4 

read4 

5 

READ5 

6 

read6 

7 

RESRVl 

8 

READS 

9 

DEPCAL 

11 

K^torj 

12 

PRINTl 

15 

DECIDE 

17 

)Z^RDER 

18 

TEST 

19 

READ? 

32 

HALFl 

53 

MATNV 

6k 

CLKRD 

72 

PRINT5 

73 

YINV 

78 

READMX 

105 

ADDMAT 

135 
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SET-UP 


JJCALC 

10 

DEPEND 

20 

MATINV 

2k 

SETUP 

27 

SIZE 

SIZERT 

iMATDIV 

k2 

110  (F(»TRAN  IV  version) 

1+4 

SETUP? 

47 

HALF2 

54 

SETUP5 

68 

SETPlU 

71 

CLKRD 

72 

PRINT5 

73 

SETPI6 

76 

YINV 

78 

SETP3 

97 

SETPll 

98 

READMX 

103 

SETPI7 

105 

SETPI5 

115 

ZGDAl 

116 

118 

SETP18 

119 

SETP19 

120 

SETP20 

121 

WCALC 

129 
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SOLVE 

All  solution  blocks  (of  programs)  contain  the  following 


JJCALC 

10 

PRIUT2 

14 

DEPEND 

20 

STEP 

21 

DEQSOL 

22 

SCCALC 

25 

MATINV 

24 

SSDET 

25 

WTSUB 

26 

STDSS 

28 

SUBl 

29 

QCSUM 

40 

MATDIV 

44 

HALF5 

6? 

CLKRD 

72 

PRINT5 

75 

YINV 

78 

PRINT? 

85 

SETSTP 

92 

PRINTS 

94 

PARTC 

95 

and  one  or  more  of  the  following  SOLV  routines : 


► 
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SOLVl 

46 

STEP 

21 

SETSTP 

92 

S0LV2 

56 

SOLV5 

57 

SOLV5 

62 

SOLVb 

48 

SOLVT 

45 

SOLVS 

77 

SOLV9 

59 

SOLVIO 

104 

SOLVll 

99 

S0LV12 

108 

SOLVI5 

109 

SOLV14 

112 

SOLVI5 

114 

ZGDA2 

117 

ZGDA3 

118 

SOLVI6 

l4l 

SOLVIT 

95 

SOLVI8 

147 

GAUS 

148 

SOME 

150 

SOMEl 

151 

ERF 

149 

SOLVI9 

161 

SOLV99 

70 

QSCALC 

152 

y 
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ITERATE 


PRINT5 

15 

PRirrr4 

16 

PRTIAL 

50 

AMX 

51 

DEVAMX 

55 

RESDET 

54 

ADORES 

55 

CEMGA 

56 

CONDET 

57 

TRMNAT 

58 

PRINT6 

Ik 

HALF5 

81 

half6 

82 

NEQS 

13k 

VIII  8 


August  1966 


FINAL  RSSULTS 

JJCALC 

ERINT2 

DEPEND 

MATINV 

CjZfRCjZl 

SIZE 

SIZPRT 

MATDIV 

HALF? 

CDCRD 

PRINT5 

YllfV 

PRINT? 

PARTC 

PUNCH 

PUNCHl 


(fortran  IV  version) 


10 

ih 

20 

2k 

kl 

k2 

no 

72 

75 

78 

85 

95 

135 

156 
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EI'ID 


DEPEND 

20 

MATINV 

24 

MAID  IV 

44 

SPCL 

55 

SPCU 

59 

SPCL3 

60 

pi/5tt 

61 

Pl/fl 

63 

SPCIA 

66 

SPCL2 

69 

CLKRD 

72 

PRIOT5 

75 

YIIW 

78 

PI^l 

84 

PI^2 

85 

PI^3 

86 

PI^IO 

87 

wiiraup 

96 

jmDPii 

100 

WNDP12 

107 

wropii4- 

122 

WNDP15 

123 

WNDP16 

124 

^^nropi? 

125 

WNDP18 

126 

V/NDP19 

127 

WTnDP20 

128 

WCALC 

129 

PUNCH 

133 

half8 

137 

INF^ 

158 

SECTION  IX 


COMMENTS  AND  SUGGESTIONS 
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COMPUTIIJG  TIME 

Problera  Solution  time  refers  to  the  total  time  necessary  to  run  a 
problem. 

Equation  Solution  time  refers  to  the  time  required  for  a single 
solution  of  the  model  equations  to  yield  values  corresponding  to  all 
T values  specified  under  "data". 

The  member  of  equation  solutions  (e)  necessary  to  solve  a problem 
cannot  be  predicted  precisely  but  will  lie  in  the  range  defined  by; 

(1+V)N  + 1 < S < (6+V)N  + 2 

where  V is  the  number  of  independently  variable  lambdas  and  N is  the 
number  of  iterations  in  the  solution.  For  zero  iterations  (simulation) 

E = 1. 

The  time  for  the  solution  of  a problem  equals  the  time  for  a 
single  equation  solution  times  the  number  of  equation  solutions  involved. 
The  time  for  a single  equation  solution  depends  on  the  model  employed. 

a)  Analytic  solution  models.  The  time  for  these 
solutions  is  usually  short  and  depends  mainly  on  the  complexity 
of  the  equations  and  on  the  number  of  data  points. 

b)  Ilumerical  solutions.  This  applies  mainly  to  the 
solution  of  differential  equations.  The  time  for  a single  equation 
solution  is  approximately  proportional  to  the  product 
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where ; LM  = maximum  value  of  X . . 

JJ 

TM  = hipest  T value  for  data 

L = number  of  lamhdas  in  prohlem 
(For  an  IBM  709^  a-  proportionality  constant  1.6  x 10  will 
yield  the  time  in  seconds.) 


IX-3 


i/ia,y  1967 


T-interrupts  may  sometimes  t>e  employed  to  reduce  computing  time 

in  the  numerical  solution  of  differential  equations.  It  requires  the 

reduction  in  the  values  of  the  largest  X . . in  the  model,  and  can  only 

J J 

he  done  if  after  such  a change  the  solution  remains  effectively  unaltered. 
Consider  the  follo\'ring  as  an  example  : 


The  equation  solution  time  of  this  model  is  governed  hy  the  largest  X . . 

J J 

which  is  Xii  = 50*  + 10 • = However,  it  is  obvious  that  compartment  1 

is  essentially  depleted  after  .1  uiits  of  time.  Thus,  a T-interimipt 
change  at  T = .Insetting  Xg^  = = 0 vrauld  not  appreciably  alter  the 

solution  and  establish  a new  X..  maximum,  namely,  X ^ = .1  + .5  = .6. 

From  this  time  on  the  solution  would  proceed  about  100  times  faster. 
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CONVERGENCE,  UNIQUMESS  AND  CONSISTENCY  ’ 

To  judge  whether  a solution  has  converged  emd,  if  so,  whether  it  | 
is  consistent  and  unique,  the  following  meas\ares  and  hints  are  useful.  j 

a)  Sums  of  squares  (SS)  for  each  "compartment"  and  total  SS  v 

i 

(including  statistical  information)  are  printed  at  the  end  of  each  j 

iteration.  It  is  necessary  for  convergence  that  the  improvement  in  SS  | 
for  the  last  iteration  "be  small  (a  few  percent  or  less  compared  to  | 

previous  iteration).  This  in  itself,  however,  is  not  sufficient  to 
guarantee  convergence. 

h)  In  each  iteration  a set  of  corrections  is  calculated  for  the 
independently  variable  parameters.  These  corrections  are  then  adjusted 
by  a factor  to  optimize  the  fit  of  the  data.  The  ratio  between  the  actual  I 
corrections  used  6ind  the  original  corrections  calculated  is  defined  CONAB 
auid  its  value  is  printed  out.  CONAB  values  close  to  \inity  at  the  end  of 
an  iteration  suggest  proper  convergence.  CONAB  values  much  less  than  i 

\uiity  suggest  that  the  model  may  be  ill-conditioned,  or  be  non-xmique 

t 

and  that  convergence  to  a least  squares  fit  may  be  very  slow.  ! 

c)  The  magnitude  of  the  correction  calculated  in  the  beginning  of 

i 

an  iteration  becomes  small  compared  to  the  value  of  the  parameter  in 
the  neighborhood  of  a least  squares  fit. 

d)  Estimated  Standard  Deviations  are  printed  out  with  the  final 
results.  These  are  determined  at  the  beginning  of  the  last  iteration.  i 
They  approximate  the  'M;rue"  estimated  standard  deviations  only  if  the 
problem  is  near  a least  squares  solution  at  the  beginning  of  the  last 
Iteration.  Standard  deviations  which  are  large  compared  to  parameter 
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I 

values  ccjuld  meem  non -uniqueness  of  the  solution  and  slow  convergence, 

e)  If  the  actual  fit  as  inspected  by  eye  seems  "good",  in  that 
the  scatter  seems  reuidom,  and  if  SS  is  comparable  to  a "reference"  SS, 
the  solution  obtained  is  near  a least  squares  solution;  it  may,  however, 
still  be  non-\anique, 

f)  If  the  fit  contains  systematic  deviations  but  otherwise  satis- 
fies all  other  least  squares  criteria  the  solution  is  inconsistent.  The 
estimated  standard  deviations  for  such  a solution  may  be  meamingless. 
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Time -dependent  (or  T -dependent) 

It  is  possitle  to  solve  differential  equations  (or  equivalent 
models)  with  time  dependent  provided  the  function  that  defines 

the  time  dependence  can  be  simulated  in  SAAI-4  as  a com.ponent  f.^  by  the 
use  of  differential  equations o A separate  compartmental  model  is  set 
up  to  generate  the  desired  function  in  one  of  its  compartments, 
say  f^  ; the  time  dependent  is  then  set  equal  to  f^  through  the 
"f -dependence 
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Delays 

Frequently  it  is  desired  to  introduce  a delay  within  a process. 

For  example,  iron  incorporated  in  a red  cell  is  not  released  until  the 
cell  "dies”.  Red  cells  have  a life  spem  of  about  120  days  and,  thus, 
the  iron  does  not  reappear  for  a time  interval  which  depends  on  the  age 
of  the  cell  when  the  iron  was  incorporated.  A ’delay*  component  is  not 
available  in  the  program.  It  may,  however,  be  simulated  by  a number  of 
compartments  in  series : 

^ r — ^ — >0 — ^ ^ 


If  a delay  t is  desired,  n compartments,  each  with  an  average 
turn-over  time  -r/n  = n/x)  may  be  set  up  in  series.  All  the  X*s 

may  be  set  equal  to  each  other  by  dependence  relations.  The  larger  the 
number  of  compartments  ^ the  closer  the  series  of  compartments  approach 
a true  time  delay.  In  a rough  way,  the  nvimber  of  ccanpartments  in  the 
chain  approximates  the  ratio  of  the  delay  time  to  the  resolution  time 
of  the  data: 


ix-8 


my  1967 


SIMUIATION 

One  can  obtain  solutions  for  a model  in  the  absence  of  any 
data  by  entering  "dummy"  data.  Such  '&ata"  need  have  entries  only  xinder 
CffltPONENT  and  T - leaving  everything  else  on  the  "data"  format  blank. 
The  program  will  assign  unit  weight  for  each  entry,  set  the  number  of 
iterations  to  0 and  proceed  with  a single  solution.  The  calculated 
values  will  correspond  to  the  entered  "data". 
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VARIABLE  INITIAL  COiroiTIONS 


There  are  no  provisions  in  the  program  for  dealing  with  variable 
initial  conditions  directly.  Indirectly,  however,  this  can  be  accom- 
plished with  a "T-interrupt"  by  introducing  "initial  conditions"  over  a 
finite — but  very  short — interval  of  time  (instead  of  instantaneously). 

A source  "feeds"  initial  conditions  into  the  desired  coupartments 

through  connecting  . which  can  be  made  variable.  The  "feed"  period, 

^ J 

is  terminated  by  setting  the  connecting  . to  zero  at  the  "T -interrupt ". 

^ J 

The  following  is  an  example 


Compartment  X is  the  source  (a  function  generator  that  is  constant) 

for  the  "feed"  and  X.  and  are  the  "feed"  parameters.  If  the 

IX  jx 

initial  condition  in  compartment  X is  f (O)  and  X.  and  X.  are  "ON" 

X XX  JX 

for  a short  interval  of  time  T,  the  "initial  conditions"  in  canpartments 
i and  j respectively  will  be  x^  f (0)t  and  X.  f (0)t.  It  is  important 
that  T be  small  compared  to  the  reciprocal  rate  constants  in  the  model. 
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PARALLEL  PATHWAYS 

Sometimes  one  desires  to  introduce  2 or  more  parallel  i)athways 
between  2 compartments,  (for  example,  the  simultaneous  presence  of 
passive  diffusion  and  active  transport).  The  program  does  not  pemit 
duplicate  parameter  designations.  This  can  be  accomplished,  however, 
indirectly  as  follows : 


to  compartment  x — a compartment  introduced  specially  for  this  purpose. 
This  ensures  that  the  right  amount  of  material  leaves  compartment  i. 


duced,  and  through  the  use  of  dependence  euid  ”f -dependence ” is  set 

equal  to  X . f . , In  addition,  the  initial  condition  f (O)  = 1 is  intro- 

duced  into  compartment  x.  It  will  be  observed  that  this  scheme  maintains 

f constant:  f (t)  = 1,  so  that  f.f  = \ .f. , 

X x'  ' ' jx  i X xi  i 

Comp  X can  serve  simultaneously  for  a number  of  extra  pathways 
between  various  compartments.  In  fact,  by  the  use  of  proper  dependence 
relations  it  can  also  serve  as  a source  of  constant  inputs  to  other 
compartments , 

Additional  pathways  from  ccaap  i "to  ^ require  independent  bypass 
compartments , 
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ZEROETH  ORDER  KINETICS 

Frequently  the  rate  of  flow  from  one  compartment  to  another  is 

constant  and  independent  of  the  eimount  of  material  in  the  compartment. 

One  way  to  accomplish  this  is  to  set  the  corresponding  . proportional 

^ J 

to  f j ^ — a special  case  of  the  f-dependence  feature : 


. 

X . f . = 


ij  J 


f . 

J 


f . = X . . 

J ij 
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The  follovj-ing  are  examples  of  various  problems  run  on  SAAM. 

Features  of  problems  and  the  SAAI'l  prograun  are  presented,  usually  in 
an  intermixed  way.  A listing  of  the  data  deck  required  for  each 
example  is  also  given.  The  first  and  last  3 rovrs  listed 

on  top  and  bottom  of  the  data  deck  are  given  for  the  lining  up  and 
recognition  of  columns.  These  cards  are  not  part  of  a data  deck. 

The  first  card  of  each  data  deck  contains  the  problem  Identification 
wMch  consists  of  three  initials  followed  by  a number  (MAI^I  002.01).  The 
integer  part  of  the  number  was  chosen  to  correspond  to  the  model  code 
for  the  problem  so  that  the  reader  may  readily  recognize  the  type  of 
of  problem. _ 
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Sum  of  Exponentials --Fit 

Fit  a sum  of  tliree  exponentials 

-''■'^t  “O'gt 

q = A^e  + A^e  + A^e  ^ 

to  the  following  data,  given  initial  estimates  = 1, 


.1,  = .0007 


Data : 


t 

.5 

1. 

lo5 

3. 

; to 

8. 

.10 

.12 

observed 

.22 

.049  i 

.0112 

0OOO33 

; .00030 

.00043 

■ .00049 

i .00054 

14. 

lo. 

lo. 

1 20. 

30. 

50. 

: .000575 

, .OOOol 

.000u3 

I.OOO05 

: 0 000675 

_ .OOOol 

Schematic : 


f 


1 


.^Ol"^ 


= e 


-Xo2'^ 


= e 


"^05^ 


% - "^l^l^l  + ^42^2  ^43^3 


Note ; The  upper  and  lower  limits  on  ^k-2’  *^43  ^01^  ^02^  ^03 

► are  somewhat  arbitrary. 

The  data  deck  is  shown  on  the  follovring  page. 


1 2 

12345u789012345o76J012345G7! 


3 4 5 b 7 8 

0ul2345b78'.3C12345b789C12345o7b'301234567S9012345G7S90 


2 saan23  nanU02.02  exanple-sun  of  exponentials  jul  Idbo 


4 3 

.Cl 

T 

X 


ICO . 


4. 

.3 

.22 

4. 

1.0 

.049 

4. 

1.5 

.0112 

4. 

3.0 

.000330 

4. 

6. 

.000360 

4. 

o . 

.00u43u 

4. 

10. 

.000490 

4. 

Xx.  • 

.000540 

4. 

14. 

.000575 

‘■r  • 

lb. 

.000610 

4. 

lu  . 

.000630 

4. 

20. 

.000650 

4. 

30. 

.000675 

4. 

50. 

.000610 

20 

X 

1. 

2 

1. 

3 

1. 

2G 

1 3. 

2 .1 

3 .0067 

2G 

4 

1 

-100. 

4 

2 

-100. 

4 

3 

-100. 

2G 

26 

2G 

. 9o  ,9b  i 


1 

n 

5. 

5. 

5. 

3 

10  G . 

100. 

100. 

4 

5 

6 


12345078y012345G7o9012343b789012345b78901234307890123450789012343G78yul234507890 

12345G78 


4 


o 
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LXAWrLt 


t>  [ n(  liuv  I 1 1 'l- 1 Ou^  • 4 


0^  i-XPONEhll/uS  JuL  lV6b 


“.uuij 

-.ouDou 

iuu.uuOiJlJ 
• uilUUb 
‘t  . uu  JUO 
4 . UDouU 
•i  . .^uuUu 

■I  »UL  JUU 
't  « kM'U  JU 
H«Ol/UUU 
•t.uUUUU 

'• . ouuuo 

•> . JOUuU 
•4  . U l<  U u u 
4 • 0(j0  Lt  J 
4 . uOOUU 
H • uOUUU 


-.i>4  00o  .tinofj 


• ji'iiOuGu 
l.cLOUUbO 
1 . jOUbUliCI 
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h.LbllUUUlJ 
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iA.uOOUUOl 
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1U<  oiiCuuUii 
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-.bunuono 
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-.IbUuUOO 
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- » dOdOuOO 
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o.uUunuOO 

-.onuouou 

iCd.OOuOoOO 
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.VrtUOb 
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“.OUOUOoO 
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-.OuObOOO 
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-.OuOOOuO 

“.OuOuOOO 
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-.OuOuooO 
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-U  -.uOUOuO 

-0  -.uOOOOO 

-U  -.UOoOUO 
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-0  -.UOliOuO 

-U  -.OUuOOO 

“0  -.oUUOuO 

-u  -.UUOOOO 


“ . 0 Cl  o 0 0 0 0 


•s\ AM\- Hi,  i^LO  fl.GL.LfV 


Mh  1 1 I I L ..I-  i<L5f  mI  Cm  at<A;  iLH 

jAAT*. 

NUiMch  Of-  CGr..POiJf  i.Tb  i*:.  4 

i^UMi.c.f-'  Of*  'jAfik  HCii'Tb  li  1h 


u/V  I A. 


.L'OouOu-o  I. 
. lUu’JOo  + ul 
. I'ibObU  + ol 

.v'UUUOC»  + Ul 
• rOJtiubO  + ol 
..•.OLiiOO  + ul 
. lOUilOU  + Uf 
. l2Uuao-*y<i 
. 14bOOU  + u^ 
. ioUUUU  + Otf 
, i4UUUl)  + bc 
.,:OoOUO  + Uw 
..500i)lUi  + 0d 
.t>tJL»OUU  + U<i 


,«:<;0U0ll-i(0 
.4SOU0P-C1 
. 1 ItuMfO-ni 
.3«.uob0-0.'^ 
. iuUoon-uJ 
.4juoor-o3 
.4VOU00-U3 
.!j40o00-03 
.b7Lo0C“03 
.oi00UC-Ci3 
.o.iOdUCi-03 
.Obl;uOii-(i3 
.r.71Su00-U3 
.oll.’J00-03 


ii.i  1 i'>«_  Lji^L’I  I lOlib 


i F(l.fi)  Jl'tyYd) 

1 .IDuOt'-*-')!  1 

£.  .lOUiJU  + oi  X 

.iOuUU  + ul 
4 .uOuOu  L 


vUl  UU) 

-♦vliObCJ  -.bCuUo 

-.(luOyO  -.OOoClU 

-.(fUObO  -.OOoUU 

.^oObO  «OUoOb 


PArA.^IlI  CKb 


AL  JUbT  AfJi»E 

liailAL  LbTlf.ATL  ‘IN 


LaI'.L.Ja(  0.  i)  = 
Lh^oUAC  U»  D- 
LAr.LUAC  U»  o)- 
( H>  1)  = 
bly^.A  I 4» 
b 1 bi'iM  ( 4 » .S ) — 


.3bU0U0(;  + Ul 
.loouuuu-ton 
.b  /U0ijnu-U2 
-.OUUOdOb 
-.oououuu 
OUUOUOo 


-.booc 

“.UOOO 
-.OuOO 
-.iuoo+ri3 
-.Ib00+u3 
1O0O  + C13 


NO  hIPlNuLLCL  rCLmTIoNShII  b 


NO  li.UcPLHULNr  STATIbllCAL  C0l,S7h>»lhT? 


. l<;H40i“03 
.24F.ott-U2 
• 2b3«i2o+0 1 
.2j7-:f9o  + 01 
. io6ol2+01 
.12fi4f.i  + l'l 
.10577J+01 
.9b2oP4-nO 
. tJLfW0.3-0n 
.777ilu-on 
.730023-00 
.67^949-00 
.620903-00 


Al  .OOOo 
.OuOuOu 
.OoOuOO 
.obOcPo 
.OOOUOU 


.5u0u+ul 
.500U+01 
.5000+ul 
.1000+U3 
. lOOU  + 03 
.luOU+03 


UP  AT  .0000 
.OuTGOO 
.OonoOO 

. OuOoOu 

.ououoo 


l/LPENOLN  I 


LiP  AT 


noouoo 

noouou 

ooouoo 

noouoo 


OP  AT  .0000 
.000000 
.000000 
.OOOOUO 
.000000 


FIXED 

VALUt 


solution 


MULLL  CODt=  2 

LbTIMAfE  OF  bIG  FKUM  P£AU-Iim  uATA  = 


77iOti71-U9 


X-.r.ev^ 


LAl-;uuA(  Ir  1)  = 
LAhDuA ( 2 » - 

Lrth.buAC  i)  = 


-.31HJUU0U+U1 
- . 1 U 1 0 C 0 U + J 0 
“.A7ij0U0U“02 


PAkA.^LIEK  k/ALUt.b 


AubUijfMbLt  l-AKAMbUI?b 


si&MA  { 
blbi'lA  t 
;j1Gi^A  ( 
LM-iuUm  ( 
LA^bUA ( 
cAtnuAl 


-b<HfaOHo9“OU 
.b72t»61H-nj 
• b7Utt9b5“0 j 
.bOOOCUOfUl 
. lliO(ll)UU  + (JU 
. u7u(J0uu-0<i 


L 

1 

2 


y 

lb 

11 

12 

li 

14 


4 


.t:0ul>ULl-oi: 
. 1 CuUbU-t  ul 
, iDUO  Jl.l  ♦01 
.AUoUOO+Ol 
.t  nuOOU  + 01 
.(  L)liOU04u1 

.UluOUll40r‘ 

.12UUUU+02 
. muuuu+Oi 

.lbUU0tJ4U2 
.180000+02 
.<•  0ij00il  + o2 

. jooo  no+u^' 

, ‘.jllOi)0U+  Ot 


JEPLUUCl.T  PAPAMtTcR! 


^C/K 

,220oby-00 
.4y2'jyp-(’l 
.lloy7i-c,l 
,3t9j»'7-03 
,3bbO00-G3 
.4o3o0'^“03 
.4y3o40-03 
.bH0'sb“-U3 
,b77o74-03 
.bUbO02-U3 
.6c7o03-03 
. hHbubr--o3 
.ticpyt'i-n3 
,M7il  1-03 


•luOuGu+nl 

.ioooOu+01 
, locuoo+ni 
. luounn  + ni 
.loOjou+ci 
. luPUOu+Ol 
. luOoOu+ui 
.luOuOu+Ol 
.lOOoOu  + C'l 
. lOOUCJu  + Ul 
■lonuOu+oi 
.lOOuOu+Ol 
. loCoCo+01 
.lOl'oOu  + Ol 


bUis  buuAPt.b  /»FTEP  U 17l‘  ATIOIjS  lb  •tioJ'^^luM— Oy 

oolip  4=  ,t:io3yii-ny 

biw  AFl£h  U ITlI  AIIONS  Ib  , 1 yi?64t)2-i ^ 


•juLUliwM  I1IL=  .lUi'JUtl 


.2t.nofly-oo 
.4'?3o9b“Ul 
.1x0971-01 
. 32Qbr»7-U3 
•bbAoPu-oS 
. 4 j3ft0b-03 
.493o4o-U3 
,5+n9by-o3 
,b/7o?H-03 
. jO^:  — 03 
, b27bOb-U3 
. bT3o^u-u3 
,6«8y0l-03 
.Ai7i  U-u3 


FIXED  I'AMAKETEHS 


00 

P20o0(J-0n 

490000-01 

112U00-01 
^jOOOO-03 
3b0000-03 
430000-03 
U90000-03 
,b40u00-03 
,875000-03 
,810000-03 
,f.30U00-03 
, AbOUOU-03 
.875000-03 
,r, 10000-03 


<jc-uc  \.c/oo 

-.Aflfl942-03  1.0031 

-.30978b-03  1.0063 

.102660-03  .9906 

.432934-06  .9967 

.200004-05  ,9944 

-.360453-05  1.0064 

-.364630-05  1.0074 

-.956971-06  1.0018 

-.287446-05  1.0050 

.369825-05  ,9939 

.219666-05  .9965 

.634175-05  .9902 

.609887-05  .9910 

-.711081-05  1.0117 


I I ti\A  I I OM  llO^'''FK  1 


L-»_  I LiNi-.il. «i.T  UF  A-i-.ATiaX=  l*ANilLOo  -.lt.4D7/^v^ 


it.r  C.<’  'A  r 1 Ol  I LOnls’t. 
i ul^lO.KLb 
1 -,b-3i..o-02 

fc  -.89154-02 

w.  .iolc./-02 


Ifb  hill  I oAlCULmIIOP  Ahu 
ui'.RtS  0(<lu.  C|. 

-.4  6 360-0  2 -.<ri  579-0  b 

— »Ho3ly— 0?  —.12843— Oo 

.15726-0?  .5..767-0? 


lOK  OF  HEj 

Lirrur  Ci<  co..-u.i-.o.  littlc.  a(i*d 

-.479a2-05  ,loOuO+Ul  .44963-03 

.lS4?l-nS  .Iu9i9+Ul  .c'3b54-o2 


.9ly94-05  ,2o073+Ul  ,Jb79b-02 


POO. VAR, ( I 
.H237R-0D 
.40471-Ou 
.64253-07 


1 ) 


LN^f.i  LCIICub  rOl<  Ai_wl>b  f /'iiLk_  »_A>-.HuAb 
(tc.b(  :)  = -.4.124  042-02 


rti_b  ( 21  — — » it? 

31  — • 15  7c.  OJ4  — 11? 


. i OoOoOtM  U 1 


P.NuAiiL  IcK  xALbcS 


AuOUb 

bibiiA  ( 4 
..>iOi-lA  ( 4 

LAhllUA(  u 
LMf.JtAf  0 


PAlvAfiL  I ullb 
I il=  ,9Ns323bb-0u 

. c)=  -.932h34c-03 

f 3)=  .94u7C55-03 

f l)=  .29951c4+0i 

f £l=  ,913o014-0x 

► 3)=  .b272M34-02 


.1tl‘c.huL..T  r>ARAMtTc.nb 


be.  bi»uAl\i_b  mF'TLP  1 ITlAmTICI.S  Ib  .10073593— Ob 
OOl-.r  4=  .1037359-09 

bxNi  Mplth  1 lTi_!</>  riOllb  lb  .94305391-11 


bui . bbuAht-b  aFTlP  1 ITLPATiOlJb  lb  .10374646—09 
COMH  4=  .iO374ub-09 

bic  mFILK  1 ITuRAlIOIlS  lb  .94314980-11 

mFTlK  a TFiLS  ll-.'  Th,L  iTn  llcUATIOlvi  - .1093?13+lU 


Soi-I  bviuAUcb  mFTLK  1 ITlFATIOIjS  lb  . 1 L)j592b4-09 
^.OHF  4=  . il)359t.5-09 

bio  AFIcU  1 ^Tl_I^A^10Nb  Ib  .94175035-11 


TuIAL  Rtb  f OR  iTEt^ATlOfj  i 
i\i_bl  1)  - -.52boM20-02 

fvubt  2)  = -.943o4602-02 

J<Lb(  31  = .17194430-02 


l(LK/\7lOi«  I IKL—  .uoOoO 


KIaFD  PARAI-ETFKS 


nLi<AI  lON  hUMhEK  2 


u(_l  cUiUHAU  I OF  A”i*iaTkIX-  A + ANIIlOo  - . 1(j5o40+ j2 


lur  OkMaTION  bOI-^ML  ol  ED  V.  ITli 
1 OKiG.KtS  i.iCn.KES 


CALCULAliOr;  AND  MCjIFKAI  10')  OF  KEb 
ORIo.  CF<  LIITLE  Ck  COkU.i.O. 


LlTTLt  A(l»l)  FOD.VAP. (I»I) 


-.51317-03 

•4ulb4-u3 

.bor>5o-04 


-.439bo-05 

.51b4o-03 

.45517-05 


,21332-09 

.52718-Ob 

.27440-07 


.47051-06 
. ift36e:-05 
..i«U65-06 


, luOOO  + Oi 
,lu9o4+01 
,3u4fl5+01 


.46338-03 

,«:5b32-02 

.47793-02 


.79005-07 

.11777-00 

.34264-00 


CoKKEC  f lOhS  FOk  AuJUbTABLL  LAi*IHuAb 
Kcb(  1)  = -.4396610-05 

KtSl  2)  = .5104794-03 

Kcb(  3)  = .4551670-05 


CUNAO 


lOOOuOO+01 


PAKAMd.r£K  VALU6.S 

AUJUSTABLt  PAHAME1ERS  DEPLNUEuT  PARAMETERS  FIXED  PARAMETERS 


SIGMA  ( 4» 

i)  = 

.9828683-00 

SIGMA  ( 

2)  = 

-.9336879-0^ 

SIGMA  C 4» 

3)  = 

.9460415-03 

lAMBUm(  U» 

1)  = 

.2994709+01 

lambda ( U » 

c)  = 

.9103002-01 

LAMdUA ( u » 

3)  = 

.6423994-02 

5Ui.i  SQUARlS 

after 

2 ITt-KATlONS 

IS  .10019912-09 

COMP 

4= 

.1001991-09 

Sib  AFTER 

ITl 

RATIONS  IS 

91090112-11 

COi'^Ad  aFTlR 

1 TRIES  l\»  The  tf-TU 

ITERATION  = .1234991+Oi 

SUM  SQUARES 

AFTER 

. d ITLKATIONS 

IS  .100070b4-09 

COMP 

4=  .i0007U«-09 

Sib 

AFTER 

2 ITlKATIONS  IS 

.90973494-11 

TOTAL  RES 

FUR  ITEKAT10i>i  l 

RlS( 

1)  = 

-.54240227-05 

KlS( 

2)  = 

.63764700-03 

RES( 

S)  = 

.56212302-05 

ITEKATiON 

T I WE-  . U lbo8 

iiekation  number  3 


ulTEKMINANI  uF  A-('iAThIX=  l*ANriLOO  -.165U0b+u2 


K^IFORHATION  CONNElTED  U17H  calculation  and  I'iOuIFICATION  of  HEb 


ORIG.RES 

MOD. RES 

Ok  I G . CR 

little  Ck 

CONQ.NO. 

LITTLl  A(I»I) 

MOD. VAR. (Irl) 

1 

-.5205b-03 

-.44446-04 

-.12114-09 

-.26729-06 

,10000+01 

.45322-03 

.35410-06 

c 

-.21995-03 

. 16O60-06 

-.24455-11 

.11535-06 

.10906+01 

.25674-02 

.86581-09 

J 

.51772-04 

.30141-05 

.lol55-0b 

.24997-06 

,3u226+01 

.48283-02 

,25243-06 

COkKECTIONS  for  adjustable  LAMBuAS 
RES(  1)  = -.‘+444810-04 

KLSI  2)  = .16o60U9-06 


KES(  J)  = .30i4lUb-05 


CUi>lAU 


.iOOOuOO+01 


PMhAMLfEK  VALUES 

rtuJUbTABLt  F-mRAMETlRS  DEPLNuEi+T  PARAMETERS  FIXED  PARAMETERS 

SiUMA  ( i)=  .Vb28b74-0U 


SIGMA  ( 

H t 

-.9367069-OS 

SiGl-iA  ( 

H f 

3)  = 

.9456869-03 

laMBGa ( 

0 » 

i)  = 

.299^663+01 

LmMBUA ( 

U t 

2)  = 

.9120155-01 

lambda ( 

u » 

0 ) = 

.6428076-02 

SUl'.  SUUARlS 
COMP 

Sib  AhlEK 
cOisAo  AFTek 


AFTEfi  3 1Tl[<aT10I^S  IS  .10uCj379-0y 
H=  .iUU053«-0'9 

3 ITLKATI0l>i5  IS  .90957990-11 

1 THlLS  IN  ThL  3Th  ITERATION  = .5G740o9+0i 


SUI.  SLUAr(uS  aF  ILR  3 ITLRATiOiJS  IS  .1000.^329-09 


lomh 

4=  .1000233-09 

SlG  AFTER 

3 itlrations  is 

.90930264-11 

TbiAL  kES 

for  ITEfvATlON  3 

KlS(  i)  = 

-.22551417-03 

RlS(  2J  = 

.64470958-06 

KcS(  3)  = 

. 15293830-04 

IIeRhTION  T1ME= 


UUCOO 


-.Ojf 

0^  L>r*ONf I’-TlALr-  JuL  ly6o 


VkloLS 


AsjLi_  h 
( Hf 
•jlbi'i/i  ( 4* 
biOr-.A  t HI 

Lm^duA(  Of 
L.Ml'l,>ly^  ( U t 


MiiAt-iL  I ( I'S 

I) -  .'rfA£7(;ii,-no 

*:)=  -.'?J'J2(JV4-0C. 

J) =  .VHouM /^-G^ 

i)=  + 

dJ=  .‘n?u2*-^-ni 

J)-  .c‘t4u3bt.-0«. 


jrPc.NL.LuT  PAPAMLTc.fb 


FIXFD  PAfcANFTFR5 


b 

7 


10 

11 

!«. 


.bC'UUGlj-uO 
.10i^»Ou-tul 
.IbUUGU+Ul 
. JOUUOU40i 
• t.OoUbo  + ul 
.r^OuGOO  + ul 
, iOu(iUU-»  0<_ 
.i2buOU+UJ 
. 14UUUU4U<; 
.16UUU04-U2 

. 1 ^ouoc+o<. 
.2nuu004o<^ 
. JULiOUu-iUf 
.DUOUOb-tOr. 


i.C/K 
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function  - Simulate  (1.02) 


Simulate  the  sine  function 
q = 1.2  sin(3t  + .5) 


Schematic 
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ii(0)  = A sin  o = 1.2  sin  .5 


— ^ = -9  fi  ^’2(0)  = A'.!)  cos  = 3*6  cos  .5 
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q = fx  = 1.2  sin  (3t  + .5) 
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Differential  Equations  - Simulate : 


q]_  = 10  f^ 


= f^  + o5  fg 

where  f^  and  f^  are  the  solutions  of 
df 

— = -5  f + 2 f + 1.2  sin(5t  + .5) 

dt 

^1^2  -2t 

— = 1 f - 2.5  fp  + 0.5  e 

dt 

f^(0)  = 0 f^CO)  = 0 

Simulate  above  for  5 units  of  time,  then  cut  off  input  function  to  f^^ 
[1.2  sin(3t  + .5)1«  Continue  solution  for  3 more  \inits  of  time  and  cut 
off  input  function  to  f^  [.5  e and  observe  solution  for  3 more 

units  of  time  after  that.  Obtain  values  for  and  q^  in  intervals  of 
.1  units  of  time  over  entire  simulation  period  and  values  for  q^  in 


intervals  of  .1  units  of  time  over  last  period  only 
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Schematic : 


G,(t)  = f^(t)  = Git)  = f^(t)  = 1 


Component  6 is  a "dummy"  introduced  for  convenience  only  as  a 
for  the  compensating  pathways. 

Gx(t)  and  G2(t)  are  function  generators. 
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Rote:  cards  with  "200."  entry  under  COMP  are  "data  generation  control  cards". 

In  this  problem  they  generate  ^ points  in  increments  of  0,1  units  of  t, 
starting  with  value  of  t just  preceding  the  control  card  (or  last  t 
value  generated  by  previous  data  generation  control  card). 
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IOWA  CITY,  IOWA  52240 


College  of  Engineering 
Department  of  Chemical  Engineering 
Area  319:  353-4875 


December  22,  1966 


Dr.  Mones  Berman 
National  Institute  of  Arthritis 
and  Metabolic  Diseases 
Bethesda  14,  Maryland 

Dear  Dr.  Berman: 


Here  is  a description  of  the  way  we 


have  adapted  your  SAAM  22  program  to  our  7044 
computer.  You  asked  for  this  in  our  telephone 
conversation  of  a few  days  ago. 


w 


Dept,  of  Chemical  Engineering 


Computer  Center 
University  of  Iowa 
Iowa  City,  Iowa 

SAAM  22  has  been  up-dated  to  Fortran  IV  language  with 
one  MAP  subroutine  for  use  on  a 7044  system.  The  number  of 
data  points  was  reduced  to  100  and  all  items  dimensioned  at 
61  were  reduced  to  50.  DS  and  CAT  were  eliminated.  These 
changes  were  necessary  in  order  to  have  the  program  fit  the 
computer.  Link  4 is  the  largest  link  and  leaves  approximately 
1090  decimal  of  unused  core. 

At  the  University  of  Iowa  we  use  an  object  deck  tape  to 
up-date  the  programs,  and  when  also  used  for  execution,  load 
time  is  about  7 minutes.  We  have  also  made  a reload  tape 
which  reduces  load  time  to  almost  zero.  The  University  of 
Iowa  will  copy  these  tapes  if  tapes  are  furnished.  To  up- 
date the  object  deck  tape  the  Delete  card  is  used. 

The  sustem  including  10  CS  uses  8 files  - FTCOO,  FTCOl, 
S.FBjZJU,  S.FBPP,  S.FBIN,  FTC03,  FTC04  and  FTCOl.  The 
following  subroutines  are  also  on  the  system: 


CHNRTN 

F04 

CNSTNT 

F05 

P0STX 

F06 

I^S 

F07 

RWD 

FPR 

ICV 

RWB 

INTJ 

ACV 

RWT 

ECV 

UTV 

HCV 

FPT 

XCV 

XBM 

FFC 

XIT 

INVERT* 

FOO 

FOl 

F02 

F03 

This  is  a U.  of  I.  sub- 
routine and  can  be  furnished. 


2 


The  Main  Line  contains  SR43  and  SCH. 

The  43  is  a NIH  number  and  SCH  is  a U.  of  I.  subroutine. 
CLKRD  is  a MAP  subroutine  from  U.  of  I. 

The  rest  of  the  links  are  arranged  as  described  in  SAAM  22 


literature  and  the  subroutine  numbers  and  names  are  kept  the 
same.  With  the  reload  tape  the  following  control  cards  are 
used : 


$ IBSYS 
$ JOB 

$ COMMENT  USE  RELOAD  TAPE  283 
$ *READTSCRATCH  TAPES  ON  B-S  AND  B-6. 
$ PAUSE  MOUNT  TAPE  283  ON  B4 
$ TIME  15 

$ PAGE  300 

$ IBJOB  NIH 
$ RELOAD  U09 


5 10  15  20  25  30  35  40  45  50 

6 7 3 2 i 4 8 8 8 


These  are  tape  assignments. 
NTAPEO  = 5,  NTAPEl  = 6,  etc. 


Data  Decks  as  used  for  SAAM  22,  9B20 , etc. 

9 16 

7 9 SAAJ4  FINAL  CARDS 

$ IBSYS 

$ CLOSE  S.SU09, REMOVE 

For  the  object  Deck  Tapes  the  following  controls  have  been  used, 

$ IBSYS 
$ JjdB 
$ TIME 
$ COMMENT 

$ PAUSE  M)Z5uNT  tape  PYB20  AS  INPUT  0N  UOS 
$ EXECUTE  UPDATE 
$ RUN  UPDATE 

$ 0BTAIN  PRINT  SUMMARY 
$ NUMBER  ZOOOOOOO 
^ $ HELETE 

W $ IB  JOB  NIH  MAP,  N)Z5  SOURCE 

$ CHAIN  SAAM  22 
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$ FILE  'FTCOO.',  UOO,  UOO,  BL0CK=15  , SINGLE,  REEL  , 
SCRTCH,  TYPE3 

$ ETC  LRL=14,  RCT=1,  E0R=  REjZJRX.  , E)Z5F=RE0FX. 

ERR=RERRX. 

$ FILE  ’FTCOl.,  U05,  U05  , BL^K=15,  SINGLE, 

REEL,  SCRTCH,  TYPE3 

$ ETC  LRL=14,  

$ FILE  ’S.FB0U’,  jZ5U,  , BL0CK=1OO  , MIXED,  TYPE3 

$ FILE  ’S.FBPP*,  PP,  PP,  BL0CK=1OO,  MIXED, SINGLE 

TYPE3  LRL=14 

$ FILE  'S.FBIN'’,  IN,  IN,  BLOCK=150 , SINGLE, 

MIXED,  TYPE3 

$ FILE  ’FTC03’,  U06 , U06 , BL0CK=15 , SINGLE, 

REEL,  SCRTCH,  TYPE3 
$ ETC  LRL=14,  .... 

$ FILE  ‘FTC041',  U07  , U07  , BL^K=15  , .... 

$ ETC  LRL=14,.  . . . 

OBJECT  DECKS 

$ ENTRY  SR89 

$ ENDCH 

5673214888 
DATA  DECKS 


9 

7 9 SAAM 


16 

FINAL  CARDS 


$ IBSYS 

$ CLjZ^SE  S.SU06 
$ IBSYS 

$ SWITCH  S.S.INl,  S.SU06 
$ STOP 
$ STOP 


For  a listing,  write  Mike  Hensel,  University  Computer 
Center,  Iowa  City.  Probably  3 tapes  would  be  required. 


« U.S.  GOVERNMENT  PRINTING  OFFICE  : 1967  0-274-091 


nih  library 


4 0125  5699 


